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insufficient  to  describe  shapes  represented  by  saooth  curves.  A  new  method  ] 
of  shape  analysis  is  described  which  exhibits  superior  performance  over 
the  critical  point  detection  methods  based  on  curvature  alone.  The  critical 
points  determined  by  this  method  are  based  on  a  set  of  coordinate  axes  that 
are  dependent  on  the  shape  itself.  This  guarantees  that  the  critical  points 
detected  are  independent  of  size,  rotation,  and  displacement  of  the  shape. 

The  results  of  applying  this  new  procedure  to  actual  shapes  are  demonstrated 
and  discussed. 

The  vector  concept  of  shape  space  is  introduced.  This  space  is  described 
in  terms  of  it's  properties.  Two  theorems  necessary  for  the  machine  recogni¬ 
tion  of  partial  shapes  are  stated  and  proved  using  shape  space  properties. 

The  critical  points  are  organized  into  structural  units  called  feature 
vectors  or  subshape  vectors  using  the  concept  of  Line  of  Sight  of  a  Point. 

The  feature  vectors  are  concatenated  to  form  a  global  shape  vector.  Shapes 
are  compared  feature  by  feature  using  a  syntactic  technique  which  will  point 
out  if  the  two  shapes  are  similar  or  not.  Examples  are  given  for  actual 
shape  data. 

Also  presented  is  a  method  for  motion  compensated  image  coding  based  upon 
a  two  step  displacement  estimation  procedure.  The  first  step  utilizes  a 
maximum  a  posteriori  (MAP)  estimator  to  determine  the  best  integer  displace¬ 
ment,  while  the  second  step  requires  solving  for  the  regression  coefficients 
that  supply  the  same  information  as  the  non-integer  portion  of  the  displace¬ 
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ones  which  use  only  linear  prediction.  This  method  is  derived  and  results 
are  presented  for  two  separate  forty  frame  digital  image  sequences.  A  sum 
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ABSTRACT 


A  system  for  the  machine  recognition  of  partial 
shapes  is  described.  Shape  analysis  metvods  are  reviewed 
in  context  to  the  problem  of  machine  recognition  of 
partial  shapes,  and  their  limitations. 

The  problem  of  defining  the  critical  points  for 
shapes  and  partial  shapes  with  various  degrees  of 
curvature  is  considered.  It  is  shown  that  the  critical 
points  derived  using  criteria  based  on  curvature  alone 
are  insufficient  to  describe  shapes  represented  by 
smooth  curves.  A  new  method  of  shape  analysis  is 
described  which  exhibits  superior  performance  over  the 
critical  point  detection  methods  based  on  curvature 
alone.  The  critical  points  determined  by  this  method  are 
based  on  a  set  of  coordinate  axes  that  are  dependent  on 
the  shape  itself.  This  gia^rantees  that  the  critical 
points  detected  are  independent  of  size, rotation,  and 
displacement  of  the  shape.  The  results  of  applying  this 
new  procedure  to  actual  shapes  are  demonstrated  and 

discussed.  '  / 

/ 

The  vector  concept  of  shape  space  is  introduced. 
This  space  is  described  in  terms  of  it's  properties.  Two 


theorems  necessary  for  the  machine  recognition  of 
partial  shapes  are  stated  and  proved  using  shape  space 
properties. 


The  critical  points  are  organized  into  structural 
units  called  feature  vectors  or  subshape  vectors  using 
the  concept  of  Line  of  Sight  of  a  Point.  The  feature 
vectors  are  concatenated  to  form  a  globlal  shape  vector. 
Shapes  are  compared  feature  by  feature  using  a  syntactic 
technique  which  will  point  out  if  the  two  shapes  are 
similar  or  not.  Examples  are  given  for  actual  shape 
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MACHINE  RECOGNITION  OF  PARTIAL  SHAPES 
USING 

SHAPE  VECTORS 

CHAPTER  I. 

INTRODUCTION 

THE  PROBLEM  AND  THE  SYSTEM  MODEL 

There  are  many  practical  applications  where  it  is 
necessary  to  identify  objects  of  shapes  in  successive 

scenes.  One  such  application  is  in  the  area  of  robot 

vision.  Here,  an  object  may  have  to  be  tracked  from  the 
instant  the  object  enters  the  field  of  view,  of  the 
robot,  to  the  time  the  object  leaves  it's  field  of  view. 
Another  application  is  in  the  data  compression  area  of 
communication  systems,  where  the  picture  data  has  to  be 
transmitted  over  a  channel  to  a  receiver.  The  picture 
data  usually  consists  of  sets  of  sequences  of  images  of 
scenes,  where  some  kind  of  motion- (activity)  has  taken 
place.  The  sequence  of  images  in  every  set  contains 

some  'common  information',  which  enables  the  ultimate 

user  (human  viewer)  to  identify  the  set  as  belonging  to 
the  same  scene.  Data  can  be  compressed  for  the  purpose 


of  efficient  transmission,  if  this  'common  information' 
can  be  identified  and  transmitted  along  with  the  first 
frame  in  the  sequence,  while  subsequently,  instead  of 
transmitting  complete  frames,  only  the  'uncommon 
information  or  changes  in  the  scene  are  transmitted.  A 
transmitter  which  is  capable  of  achieveing  this  goal 
will  require  a  'smart  receiver'  which  would  be  capable 
of  combining  the  common  and  uncommon  information  and  be 
able  to  simulate,  predict  or  reconstruct  the  pictures  of 
the  scene  to  the  desired  human  fidelity  criterion. 

The  above  process  of  extracting  information  can  be 
divided  into  three  steps  or  subprocesses,  which  are, 

1)  The  correspondence  process  (  the  extraction  of  the 
common  information  )  [88], 

2)  The  interpretation  process  (  the  extraction  of  the 
uncommon  information)  [88],  & 

3)  The  predicition  and  simulation  process. 

For  a  communication  system  the  first  two  are 
exclusively  the  domain  of  the  'smart  transmitter',  while 
the  last  is  that  of  the  'smart  receiver'.  In  a  robot 
vision  system,  the  functions  of  the  '3mart  transmitter' 
and  the  'smart  receiver'  may  be  accomplished  at  the  same 


The  correspondence  process  or  problem  is  that  of 
identifying  *  portion  of  changing  visual  array  as 
representing  the  same  scene  or  object  in  change  or 
motion. 

The  extraction  of  the  uncommon  information  about  the 
sequence*  or  the  interpretation  process  can  be: 

1)  a  qualitative  process  or 

2)  a  quantitative  process 

The  process  is  a  qualitative  process  when  the 
extraction  of  information  is  based  on  statistics  derived 
by  operations  on  adaptive  blocks  or  blocks  of  variable 
size,  within  the  picture  and  the  final  output  results  in 
statements  like,  '  the  object  located  near  the  left 
Conner  of  the  picture  went  through  a  rotation  of  1.57 
radians  about  the  z  axis  while  the  background  moved 
forward  by  about  a  feet*. 

The  process  is  a  quantitative  process  if  the 
extraction  of  information  is  based  on  statistics, 
derived  by  operations  on  predefined  or  fixed  blocks  or 
sub-blocks  of  the  picture,  e.g.  operations  based  on 
rows  or  columns  .  The  representation  of  such  information 
is  not  qualitative.  Usually  such  a  representation  is  in 
the  form  of  numbers  of  data  arrays  which  cannot  be 
directly  interpreted,  e.g.  the  2-D  array  representating 


the  difference  in  the  pixel  intensities  between  two 
frames  in  the  sequence. 

The  correspondence  process  begins  before  the  object 
has  completely  entered  the  field  of  view  and  continues 
even  while  the  object  is  undergoing  occlusion  or  leaving 
the  field  of  view.  This  suggests  that  a  correspondence 
process  based  on  global,  shape  or  object  recognition 
methods  alone  would  not  be  adequate  for  such  a  system. 
In  this  dissertation,  a  system  which  is  capable  of 
identifying  an  object,  when  it  is  partly  or  partially  in 
field  of  view,  with  itself  at  some  other  time,  when  it 
is  completely  in  the  field  of  view,  will  be  called  a 
Partial  Shape  Recogntion  System.  In  such  a  system 
recognition  of  shapes  is  based  upon  identifying  elements 
or  correspondence  tokens  in  different  views  as 
representing  the  same  shape  in  at  different  times,  and 
thereby  maintaining  the  perceptual  identity  of  objects 
in  motion. 

The  correspondence  tokens  which  are  matched  by  the 
correspondence  process  [88]  include  critical  points, 
edge  fragments,  bars,  small  blobs  etc.  These  tokens  are 
detected  first  and  organised  hierachally  into  more 
structured  forms,  and  finally  into  distinct  objects. 


information  necessary  to  alter  the  degree  of  the 
interpolating  curve  can  only  be  obtained  after  the  shape 
data  has  been  processed  at  least  once.  The  interpolated 
curve  is  then  examined  for  evaluation  of  critical  points 
and  curvature.  The  combination  of  these  two  factors 
helps  in  the  extraction  of  features  and  the  length  of 
segments  to  be  matched.  These  features  may  be  altered 
deleted  or  merged  based  on  the  knowledge  from  the  data 
base  or  the  feedback  from  the  output  table.  Though 
often,  features  may  be  deleted  based  on  their  relative 
size,  the  rules  for  altering  the  feature  structures  are 
not  well  established  and  will  not  be  addressed. 

The  features  are  then  converted  into  shape  vectors 
and  stored  in  the  data  table.  Information  about  the 
segments  of  required  length  is  also  stored  in  the  shape 
data  table.  This  information  can  be  in  the  form  of 
curvature  of  the  segments,  length  of  segments, 
interconnections  between  segments,  frequency  domain 
information  about  curvature,  and  or,  other  space 
invariant  properties. 

The  shape  data  table  which  also  includes  information 
about  block  processes  is  then  compared  with  the  shape 
data  table  simulated  by  the  predictor  with  knowledge 


Comparison  with  the  latter  helps  in  extraction  of 
symmetry  and  other  aspects  of  the  shape.  The  results  of 
comparsion  are  stored  in  the  output  table.  These  results 
are  fed  back  into  the  various  blocks  of  the  system  shown 
in  Fig.  1.  The  feedback  is  mainly  to  establish 
relationships  with  other  shapes  on  the  scene,  to  delete 
relatively  less  important  critical  points  or  features, 
to  update  the  data  base,  and  to  the  change  the  degree  of 
the  interpolating  curve  if  necessary. 

At  the  present  the  methods  of  recognizing  shapes 
[20],  [21],  [60],  [65],  [63],  [66],  [91],  can  be 
categorized  as  either  global  or  local  in  nature.  Within 
the  class  of  global  shape  analysis  methods,  there  are 
two  categories  that  under  certain  circumstances  possess 
the  ability  to  recognize  complete  or  whole  shapes 
independent  of  size,  rotation,  or  location.  These  are 
the  Fourier  descriptors  methods  and  the  Syntactic  or  the 
Graphical  methods.  The  Fourier  descriptors  based 
algorithm  performs  satisfactorily  on  complete  shapes. 
The  Fourier  coefficients  extracted,  are  indeed 
independent  of  size,  rotaion,  and  location  when  the 
shapes  are  complete.  However  this  method  does  not 
perform  satisfactorily  and  in  fact  fails  entirely  when 
the  class  of  shapes  is  allowed  to  include  incomplete  or 
partial  shapes.  An  example  is  presented  in  chapter  II 


that  demonstrates  that  Fourier  descriptors  method  fails 
to  work  on  incomplete  shapes.  The  results  of  the 
experiment  are  discussed  in  order  to  point  out 
specifically  why  the  algorithm  cannot  perform 
satisfactorily  on  partial  shapes. 

The  other  class  of  global  methods  namely  the 
Syntactic  methods  [21] ,  have  restricted  use  in 
recognition  of  shapes  because  these  algorithms  tactily 
assume  a  priori  that  the  shapes  have  been  identified  by 
their  parts.  These  algorithms  then  investigate  the 
relationship  between  the  various  parts  of  the  shape.  A 
human  shape,  for  instance  has  a  hand  or  a  face  at  some 
definite  orientation  and  location  with  respect  to  each 
other. 

The  Local  category  of  shape  analysis  algorithm  [91], 
uses  curvature  as  a  criterion  for  detecting  the  peaks 
and  valleys  of  a  shape.  These  peaks  and  valleys  are 
called  the  local  shape  descriptors.  This  shape 
comparision  algorithm  is  not  independent  of  rotation.  In 
Chapter  II  specific  examples  are  given  that  demonstrates 
that  the  present  algorithm  is  not  independent  of 
rotation. 

In  Chapter  III.  the  concept  of  curvature  is  presented 
form  the  point  of  view,  of  differential  geometry.  The 


concept  is  presented  so  as  to  determine ,  why  the  local 
shape  descriptors  in  Chapter  XI  are  not  independent  of 
rotation.  Further  it  is  shown  that  the  critical  points 
found  by  using  curvature  as  a  criterion  are  insufficient 
to  describe  many  shapes.  This  is  demonstrate  by  an 
example. 

In  Chapter  IV  the  concept  of  Line  of  Sight  of  a  Point 
(LSP)  and  Line  of  Sight  of  a  straight  line  Axis  (LSA) 
are  introduced.  These  concepts  are  then  used  to  define 
to  Adaptive  Line  of  Sight  method  (ALS)  for  determining 
the  critical  points.  The  effect  of  this  algorithm  on 
actual  shape  data  is  presented.  The  critical  point 
detected  by  this  method  are  very  close  to  those 
perceived  by  the  human  vision  system  in  most  cases.  The 
ALS  method  can  also  be  used  as  a  basis  for  detecting  the 
axes  of  symmetry  in  a  shape,  if  any  exist.  However,  such 
a  task  can  be  only  be  achieved  at  a  post  cognitive 
level;  i.e.  the  critical  points  found  by  the  ALS  method 
can  be  used  to  locate  the  axes  of  symmetry.  Location  of 
such  an  axes  is  an  important  part  of  shape  analysis, 
because,  it  seems  that  in  the  human  vision  system  ,the 
critical  points  located  with  respect  to  an  axis  of 
symmetry  are  considered  more  important  than  critical 
points  located  with  respect  to  any  other  axis. 
Unfortunately  the  ALS  method  cannot  always  detect  an 


axis  o£  symmetry  at  the  precognitive  level. 

In  Chapter  V  concepts  from  several  areas  of  shape 
analysis  [19],  [52],  [91],  are  combined  with  some 
entirely  new  concepts  concerning  shapes  for  the  purpose 
of  providing  the  foundation  for  a  new  approach  to  define 
shape  as  vectors  in  an  appropriate  space.  The  properties 
of  this  space  are  stated  definitively,  after  the  concept 
of  size  variable  has  been  solidified.  This  vector  space 
is  called  the  shape  space.  Two  theorems  usefuls  in  the 
partial  shape  recognition  problem  are  stated  and  proved 
utilizing  shape  space  properties. 

In  Chapter  VI  the  concept  of  Line  of  Sight  of  a  point 
is  used  to  organize  critical  points  into  feature  vectors 
in  the  shape  space.  These  feature  vectors  are  then 
concatenated  to  form  shape  vectors.  Procedures  and 
tests  necessary  for  comparing  shape  vectors  are 
examined.  The  comparison  procedure  is  based  on  a 
Syntactic  method  which  will  point  out  whether  one  shape 
is  part  of  a  more  complex  shape,  or  whether  the  shapes 
are  totally  dissimilar. 

The  conclusion  and  discussion  in  Chapter  VII  places 
into  perspective  the  overall  effectiveness  of  methods 
for  analyzing  shape  based  on  critical  points,  addresses 
the  question  of  thresholds,  the  problem  of  locating  axes 


CHAPTER 


OTHER  METHODS  AND  THERE  RELEVANCE  TO  THE  PROBLEM 
II. a  FOURIER  DESCRIPTORS  METHOD 

Theca  exist  two  types  of  Fourier  descriptors.  The 
first  type  of  descriptors#  used  by  Zahn  and  Roskies 
[98]#  have  been  called  descriptors  Sn  by  Pavlidis  [60]. 
In  the  method  of  descriptors  S^,  the  shape  is 
represented  by  the  continuous  function# 

a(  t(k)#k)M(  k  )  +  t(  k-  )  (II. 1) 

where  t(  k  )-  2*A(k)/L 

i(  k  )»  arc  length  between  the  starting 
point  and  the  k  th  point  on  the 
curve. 

♦(  k  )■  net  amount  of  angular  change 

between  the  starting  point  and 
the  k  th  point  on  the  curve. 

L  ■  the  perimeter  of  the  curve 

The  descriptors  ^  for  a  continuous  shape  are  then 
defined  as# 

2ir 

S  ■  lfa(t(k),  k)  exp(-j2irnt)dt 


(II. 2) 


These  descriptors  exhibit  some  notable  shortcomings. 
Among  these  are  1}  the  property  of  the  closure  of  the 
curve  is  not  preserved  2)  simple  shapes  such  as  squares 
and  triangles  cannot  be  distinguished  from  one  another 
when  only  the  vertices  are  given.  The  second  type  of 
Pourier  descriptors,  namely  the  descriptors  Tn  [60], 
[63],  [66]  exhibit  charactersitics  that  are  superior  to 

the  descriptors  S  in  the  sense  that  the 

reconstruction  of  the  shape  from  a  finite  set  of 
coefficients  leads  to  a  closed  curve.  Also  the 
convergence  properties  are  superior. 

For  the  descriptors  T  the  shape  data  is  represented 

n 

in  the  complex  form, 

u(M  -  x(i)  +  j  y(M  (II. 3) 

where  (x(£)  ,  y(t))  are  the  coordinates  of  the  point  on 
the  curve  and  *  is  the  arc  length  from  the  defined 
starting  point.  The  descriptors  T Q  are  then  given  by, 

^»llu(  t )  exp(-j2irrt  /L)  dl  (II. 4) 

•-  K 

The  Fourier  descriptors  based  algorithm  normalizes 
for  position  by  setting  TQ  to  zero  .  Normalization  for 
scale,  rotation,  and  starting  point  of  a  contour  is 
achieved  by  multiplying  the  n  th  coefficient  by 


s  exp(  j  (  4  +n  a ) )  .  The  parameter  s  scales  the  shape  to 
the  normalized  size,  the  parameter  (  $+  n  a)  rotates 
the  contour  to  the  normalized  position.  The 
normalization  is  such  that  the  coefficients  T+1  and 
are  pure  imaginary  numbers  and  their  sum  has  magnitude 
one. 


It  is  easy  to  see  why  such  a  normalization  is 
necessary  for  discrete  data.  In  this  case,  the  Fourier 


descriptors  are  given  by, 
N-l 


T(  n  ) 


*2 


•“IT  Si  x(  k  )  +  j  y{  k  ) ) exp(  -j2*n)c/N) 

k*°  (II. 5) 


Multiplying  (I.b.5)  by  s  exp-jC*  +  n®  )  the  normalized 


coefficients  ,  T  (  n  )  are  then  given  as. 


N 


T  (  n  )  ■  T(  n  )  s  expj  (H  n®  ) 
N 


(II. 6) 

Next  imposing  the  requirement  on  (I.b.6) ,  that  the 
coefficient  be  purely  imaginary  at  n  »  *  1  leads  to  the 
condition  that, 

2  jHc  -U  +  ®  )  ■  2m  +  {-  Z*  k  -(  *  -  a  ) )  (II. 7) 


where  m  is  an  integer. 

Setting  the  real  part  to  zero  by  appropriate  choices  of 
9  sad  ®  is  equivalent  to  normalizing  for  rotation  and 
starting  point  of  the  contour.  As  a  matter  of  fact,  the 
two  equations  could  have  been  set  equal  to  any  constant. 


The  imaginary  parts  of  the  normalized  coefficient  at 


n  ■ ±  1  are, 


N-l 


IM(T  (  1  ))«s  /(-x(k)sin(2*k  -  (♦+«)) +y(k)  cos  (2ir_k  -  (♦+»)) 

N  N 

(11.8) 


N-l 


IM(T  (-1  ))^s\(-x(k)sin(2*Jc  -(*-«) )  +  y(k)  cos(2wk_ -(  ) ) 

N  N L-t  N  N 


k-0 


(II. 9) 


Summing  and  equating  the  magnitude,  of  the  normalized 
coefficients  at*!,  to  one,  yields. 


N-l  -1 

a^j(x(k)sin(  A  +  B) -sin(A-C) )  +(y (k)  (cos (A-C)  +cos (A-B) )  n| 


k-0 


(11.10) 


where  A  ■  2ir  k/N 
B  •  ♦  +  a 
C  -  ♦  -  a 


Another  function  [64]  which  has  also  been  used  for 
normalization  is  the  standard  deviation  of  the  data, 


JL 


\lx(k)-i}2  + 


(II. 11) 


where 


It  is  apparent  that  both  s  *  and  o  are  homogeneous 
function  of  order  one  of  the  data  points. 

In  a  real  shape  recognition  application  it  is  not 
known  a  priori  whether  the  shape  under  examination  is  a 
part  of  a  more  complex  shape  or  a  shape  in  its  own 
right.  The  above  descriptors  in  their  present  form  are 
not  capable  of  recognizing  that  a  partial  shape  may  be 
part  of  a  more  complex  shape.  An  example  is  presented  to 
demonstrate  the  validity  of  this  contention.  The  shapes 
used  are  shown  in  Fig.  2.  Figure  2-a  is  a  complete 
shape ,  namely  a  swept  wing  plane ,  while  Figure  2-b  is  a 
partial  shape ,  namely  the  front  part  of  the  plane.  The 
plots  of  the  real  and  imaginary  part  of  the  normalized 
Fourier  descriptors  obtained  by  using  (11.10)  for  the 
complete  shape  are  shown  in  Fig.  3-a.  The  corresponding 
plots  for  the  partial  shape  are  shown  in  Fig.  3-b. 
Similar  data  obtained  using  (II. 11)  is  shown  in  Fig.  4. 
It  is  apparent  that  comparing  the  two  sets  of  data 
yields  nothing  more  that  a  statement  that  the  two  shapes 
are  dissimilar.  Two  explanations  for  this  are  1)  the 


FOURIER  DESCRIPTORS 


FIG.  4.  FOURIER  DESCRIPTORS  CALCULATED  US  DC  TOE 


parameters  s *  and  a  are  not  independent  of  the  shape 
and  2)  the  Fourier  descriptors  method  compares  the 
'frequencies'  of  the  shapes.  The  frequencies  of  the 
partial  shape  are  not  the  same  as  the  frequencies  of  the 
complete  shape. 

Besides  the  above  disadvantages  there  is  a  very 
important  aspect  of  the  Fourier  descriptors  that  seems 
to  have  skipped  the  attention  of  most  researchers  in  the 
field.  The  fundamental  reason  for  comparing  shapes  using 
Fourier  descriptors  lies  in  the  fact  that  the  Fourier 
coefficients  for  a  given  shape  are  unique;  however  there 
are  many  transforms  that  will  give  a  unique  set  of 
coefficient  for  a  shape.  Considering  the  computational 
power  of  the  present  day  computers  there  is  no  reason  to 
treat  other  transforms  as  second  class  citizens. 
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II. b 


LOCAL  SHAPE  DESCRIPTORS 


Curvature  has  been  used  to  determine  representative 
points  on  a  shape  by  many  researchers  [14] ,  [91]. 

Usually  curvature  is  defined  as  an  operation  on  three 
points  in  a  sequentially  ordered  list. 

The  local  category  of  shape  analysis  algorithm  [91] 
uses  'curvature'  as  a  criterion  for  detecting  the  peaks 
and  valleys  of  a  shape.  These  peaks  and  valleys  are 
called  the  local  shape  descriptors.  This  shape 
comparision  algorithm  is  not  independent  of  rotation.  A 
specific  example  is  given  that  demonstrates  that  the 
present  comparision  type  local  shape  analysis  algorithm 
is  not  independent  of  rotation.  The  definition  of 
curvature  as  given  in  [91]  is  examined  and  used  to  find 
peaks  and  valleys  (critical  points)  of  a  simple  shape 
(cardioid) ,  before  and  after  rotation. 

The  Local  Shape  descriptors  method  has  been  used 
[91]/  for  comparing  shapes  stored  in  a  library.  In  this 
method/  the  angle  function  is  defined  asr 
s(k)»arctan( (y(k)-y(k-l) )/(x(k)-x(k-l) )  (11.12) 

Curvature  is  then  defined  as  the  derivative  of  the  angle 
function/  that  is  s' (  k  ) .  The  peaks  and  valleys  in  s' (M 
are  used  for  finding  the  peaks  and  valleys  of  the 


curve.  Each  local  shape  descriptor  then  consists  of  two 
adjacent  peaks  and  a  valley  (  alternate  angles  and  a 
distance) . 

It  is  apparent  from  (11.12)  that  this  procedure 
results  in  a  nonlinear  transformation  on  the  data,  also, 
the  function  is  not  finite  at  *  1.57  radians.  It  may  be 
generally  stated  that  this  function  may  be  used  to 
locate  peaks  and  valleys  where  :  1)  the  angle  between 
any  three  adjacent  points  is  acute  and  2)  the  shape 
becomes  parallel  to  the  y-axis.  These  are  equivalent  to 
the  points  where  the  angle  becomes  ±  1.57  radians  with 
respect  to  the  x-axis.  For  example  consider,  the 
cardioid  shown  in  Fig.  5-a,  this  shape  was  generated  by 
sampling  at  constant  intervals  of  (6.28/256)  radians 
rather  than  the  arc  length.  This  prevents  an  acute  angle 
from  occuring  between  any  three  adjacent  points  on  the 
shape.  The  cardioid  was  then  rotated  by  1.57  radians. 
The  corresponding  plots  for  the  angle,  the  derivative  of 
the  angle  function,  and  the  peaks  obtained  using  this 
method  are  shown  in  Figures  5-a,  5-b,  6-a,  and  6-b, 
respectively.  Note  that  the  peaks  obtained  by  this 
method  are  not  unique  and  are  dependent  on  the  rotation. 
Every  different  rotation  will  yield  a  different  set  of 
peaks.  Spurious  peaks  will  occur  at  the  discontinuities 
of  the  arctan  function  that  may  be  *  1.57  radians  or 
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II. c 


MISCELLANEOUS  OTHER  METHODS 


Several  methods  have  been  utilized  in  the  past  for 
determining  critical  points  [14] , [19] , [321 •  Most  of 
these  methods  are  based  on  operations  on  a  fixed  number , 
of  points.  These  methods  were  originally  meant  for 
ideal  curves  (noise  free  curves) ,  where  the  operations 
were  defined  on  a  set  of  adjacent  points.  To  by  pass 
the  effect  of  noise,  and  to  arrive  at  a  better  estimate 
of  the  parameters  being  evalued,  these  operations  were 
redefined,  [14] ,  [19] ,  for  points  of  the  shape  located 
at  some  fixed  distance  (instead  of  adjacent  points) . 
Since  it  is  unlikely  that  an  intelligent  machine  will 
have  a  priori  knowledge  about  the  size  of  the  shape  to 
be  analyzed  or  the  relationship  of  the  number  of  sample 
points  to  the  size,  the  ambuiguities  involved  in 
detection  of  critical  points  using  methods  that  are 
totally  dependent  on  operations  on  fixed  number  of 
adjacent  points,  are  high.  A  brief  description  of  some 
of  these  methods  which  are  represented  here  by 
operations  on  three  adjacent  points  is  presented  next. 


Il.C.i  THE  METHOD  OP  CENTRO IDAL  VECTORS 


In  this  method  [19] ,  the  ith  point  on  a  shape  at  a 


vector  distance 

di  Zp, 

from 

a 

reference 

point 

(typically  the  shape 

centroid 

) 

is 

said  to 

be  critical 

with  respect 

to  the 

reference 

point 

if 

(( 

di-l- 

d1  ) 

and 

( 

d1 

-  d1+l 

)) 

or 

(( 

*1-1  ~ 

*  1  ) 

and 

( 

*1 

“  *1+1 

)) 

have  opposite  sign. 


This  operation  is  very  local  in  nature  and  is 
extremely  sensitive  to  noise.  Round-off  or  truncation 
errors  also  have  a  deleterious  effect  on  the  operation. 
It  also  fails  very  often  when  dealing  with  smooth  curves 
or  in  situations  where  the  centroid  is  located  away  from 
the  shape  boundary  as  shown  in  Fig.  7. 


Il.C.ii  METHOD  OP  CURVATURE  VECTORS. 


This  method  is  basically  a  two  pass  process.  In  the 
first  pass/  the  i  th  point  on  the  shape  with  a  curvature 
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CRITICAL 


FIG. 


7  .  A  CASE  AGAINST  THE  QCTOKAL  ME2HCD 


is  set  as  a  critical  point  if 


( 


*1  ’ 


'i-i 1  01  <  '< 


*1+1  * 


have  opposite  sign.  In  the  second  pass,  the  points  at  a 
maximum  distance  from  the  straight  line  joining  every 
two  successive  critical  points  is  set  as  critical.  This 
method  like  the  previous  method,  is  very  sensitive  to 
noise.  Though  it  works  for  most  of  the  smooth  curves, 
there  are  instances  where  it  fails.  For  example,  this 
method  produces  only  one  critical  point  for  the  cardioid 
shape  of  Fig.  5. 


IX.c.iii  THE  MAGNITUDE  OF  CURVATURE  METHOD 

In  the  Magnitude  of  Curvature  Method  the  ith  point  is 
determined  to  be  critical  if, 

!  i  *  Threshold 

The  critical  points  were  determined  by  using  this  method 
on  the  elephant  shape,  which  was  obtained  by  non- 
uniformly  sampling  a  hand  sketched  figure.  These 
critical  points  are  shown  in  Fig.  8.  It  is  apparent  that 
some  critical  points  were  missed.  For  example  the 
critical  points  for  the  trunk  are  missing.  This 


procedure  is  more  immune  to  the  noise  than  the  others; 


CRITICAL  FOINIS  FOR  TOE 


however  ,  an  equally  critical  problem  is  added. 
Specifically/  a  threshold  must  be  determined  a  priori. 
Additionally  the  method  fails  completely  on  smoothly 
varying  curves  such  as  the  cardioid  of  Fig.  5. 


II. C. iv  THE  LINE  OF  SIGHT  METHOD. 

A  procedure  that  perforins  well  for  polygonal  shapes  and 
curves  with  very  few  boundary  points  is  called  the  Line 
of  Sight  method.  This  method  has  not  appeared  before  in 
literature  and  is  being  presented  for  the  first  time.  It 
is  described  here  to  give  a  feeling  for  the  problem  of 
locating  critical  points  on  a  shape.  It  was  also  the 
method  which  led  to  the  developement  of  the  concepts  of 
Line  of  Sight  and  the  Adaptive  Line  of  Sight  method 
described  in  subsequent  chapters. 

In  this  method/  if  d7-  denotes  the  normal  distance  of 
the  i  th  point  from  a  straight  line  L,  then  the  i  th 
point  on  the  shape  boundary  is  said  to  be  critical  with 
respect  to  the  straight  line  L  if 

(  d.  .  -  d.  )  and  (  d.  -  d.  ) 
i-i  i  i  i+1 

have  opposite  sign.  The  set  of  critical  points  found 
with  respect  to  the  set  of  tangent  lines  drawn  at  all 


points  on  the  shape,  will  then  be  called  the  critical 
points.  It  is  obvious  that  for  shapes  other  than 
polygonal  an  infinite  number  of  axes  or  tangent  lines 
are  required.  Nonetheless,  the  attributes  exhibited  by 
this  method  are  very  desirable.  It  is  however  necessary 
to  reduce  the  dimensionality  of  the  problem;  this  is  the 
subject  that  is  addressed  by  the  Adaptive  Line  of  Sight 
method. 


CHAPTER  III. 


III. a  THE  NEED  FOR  DETECTING  CRITICAL  POINTS. 

The  classical  method  of  comparing  curves  or  shapes  in 
space  using  the  geometric  invariant  parameters, 
curvature  and  torsion  is  examined  here.  This  is  done  to 
demonstrate  a  need  for  detecting  critical  points  based 
on  criteria  other  than  curvature,  and  to  pin-point  the 
difference  between  the  definition  of  curvature  as  used 
in  the  field  of  differential  geometry  and  the  definition 
resulting  from  (11.12)  The  weakness  of  using 
curvature  as  a  criterion  [91] ,  for  comparing  practical 
shapes  is  briefly  mentioned.  The  cardioid  shape  is  used 
again,  to  show  that  even  with  the  differential  geometric 
definition  of  curvature,  curvature  cannot  be  used  as  a 
single  entity  for  detecting  critical  points  which  may  be 
used  for  comparison  of  shapes. 

III. b  CORVES  IN  THREE  DIMENSIONAL  SPACE 

The  problem  of  shape  recognition  is  analogous  to 
recognition  of  curves  in  space.  Therefore,  well  known 
concepts  and  theorems  from  differential  geometry  can  be 


utilized  in  shape  analysis.  Let  c  be  a  curve  m 
Euclidian  space  R  n  of  class  C  >■  1  whose  domain  is  I 

s 

where  s  is  the  parameter.  Then  a  differential  geometry 
theorem  that  is  particularly  useful  is  [40], 

n 

Theorem:  Every  regular  curve  c:I  R  can  be 

parameterized  by  its  arc  length. 

n 

In  other  words,  given  a  regular  curve  ctl^R  there 
is  a  change  of  variablese:^  ■*  Ig  such  that  |  (c.e  ) ' (s)  ]»1, 
where  (c.9)is  a  composite  function. 

Let  c  ■  c(  s  )  be  the  parametric  representation  of  the 
curve  under  analysis  with  s  as  the  natural  parameter (i.e 
i  dc/ds  I  •  1)  then  the  vectors  t,  n,  b,  satisfy  the 
Serret-Frenet  equations  [26],  [40], 
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where, 

*  ■  curvature, 

T  »  torsion, 
t  ■  tangent, 

n  ■  normal  at  the  point, 
b  >  binormal  at  the  point, 


and  the  dot  notation  denotes  the  deravative  with  respect 


to  the  natural  parameter  s. 


The  parameters  k  and  T  are  geometric  invariants  and 
their  existence  and  uniqueness  is  guaranteed  by  the 
following  theorem  [40] , 

Theorem  :  Let  k  {  s  )  and  T  (  s  )  be  arbitrary 
continious  functions  for  a  <  s  £  b  .  Then  there  exists, 
except  for  position  in  space,  one  and  only  one  curve  c 
for  which  k  (  s  )  is  the  curvature,  t  (  s  )  is  the 
torsion  and  s  is  the  natural  parameter. 


When  the  curve  under  analysis  lies  in  a  plane,  the 
torsion  t  is  equal  to  zero  and  the  binormal  b  is 
constant.  Thus  for  a  curve  in  the  x-y  plane  (III.l) 
reduces  to, 
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Now  if  e  is  the  angle  made  with  respect  to  the  x-axis  by 
the  tangent  to  the  curve  then, 
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(III. 3) 


where  i  and  j  are  unit  vectors  in  the  x  and  y  directions 
respectively,  differentiating  (III. 3)  with  respect  to 
the  natural  parameter  s  yields, 
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Substituting  equation 
yields  the  result, 
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(III. 3) 


in  equation 


(III. 4) 
(III. 4) 


(III. 5) 


Comparing  (III. 2)  with  (III. 5) 

e 

K  a  9 


(III. 6) 


In  view  of  the  different  representations  of  the  word 
curvature  [14],  [91]  it  is  necessary  to  emphasize  that 


the  derivative  in 

(III. 6} 

is  with  respect 

to 

the 

natural  parameter 

s  and 

not  with  respect 

to 

some 

arbitrary  distance 

measure.  If  a  curve 

is 

not 

represented  in  terms 

of  the 

natural  parameter. 

but 

some 

other  real  valued  function  (saye  *e  (  s  )  )  then  this 

transformation  should  be  allowable.  The  implication  is 

that,  0:1  -*•  I  is  a  monotonic,  injective  mapping  of  the 
S  0 

interval  I$  onto  I0  .  Where  Is  and  I0  are  the 
respective  domains  in  s  and  e  over  which  the  curve  is 
defined. 


Therefore  the  function  arctan(  9  )  with  values 
between  ±  1.57  radians  or  arccos(  e)  with  e  between  0 


and  3.14  are  not  allowable  changes  of  parameter  and  any 


property  based  on  these  transformations  may  not  be  a 
property  of  the  curve  but  a  property  of  the 

representation.  Thus  if  a  curve  or  a  shape  has  been 
represented  in  terms  of  the  sample  number  *  Ic *  and  if  the 
algorithm  is  unable  to  affect  the  one-to-one 
transformation  described  above,  then  the  following  fact 
from  the  Fundamental  Theorem  of  calculus  should  be 
exploited, 

(III. 7) 

□sing  (III. 7)  it  can  be  shown  [40]  that  the  magnitude 

of  curvature  can  be  obtained  from  the  following 

relationship, 

- 

|  *  |  «  ,  c'  x  c' '  I  /  |  c*  |3  if  c'  j*  0 

.0  otherwise  (III.  8) 

Where  the  symbol  x  denotes  the  vector  cross-product  and 
the  symbol  *  denotes  differentiation  with  respect  to  k. 
Curvature  is  a  vector  quantity  and  it  points  in  the 

direction  of  the  normal  to  the  curve. 

An  investigation  of  (III. 2)  suggests  that  curature 
alone  is  sufficient  to  uniquely  describe  curves  or 
shapes  in  a  plane.  The  problem  of  comparing  planar 
shapes,  is  then  equivalent  to  comparing  their 

curvatures.  This  solution  works  well  for  ideal  curves. 
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However, it  is  ruled  out  as  a  sole  model  to  be  used  by 
systems,  or  algorithms,  to  analyze  shapes  in  a  manner 
similar  to  the  human  vision.  Explicitly  the  following 
draw-backs  make  this  solution  unacceptable. 

1)  Curvature  as  defined  above  (  proportional  to  the 
derivative  of  the  tangent  )  is  an  operation  on  three 
adjacent  points  and  any  errors  due  to  sampling  or  noise 
tend  to  make  the  curvature  functions  of  any  two  shapes, 
otherwise  alike  for  most  visual  purposes, 
unrecognizable. 

2)  The  curvature  function  does  not  automatically  render 
itself  useful  to  syntactic  comparision  methods  which  is 
an  inherent  aspect  of  human  vision. 

3)  Curvature  is  inversely  proportional  to  the  size, 
which  implies  that  the  plots  of  the  logarithm  of  the 
curvature  of  two  identical  shapes  of  different  sizes 
differ  from  each  other  by  an  additive  constant.  This 
does  not  seem  to  be  a  problem  at  first  glance,  however, 
it  is  a  major  problem  when  comparing  partial  shapes  of 
different  sizes  and  different  number  of  samples.  Since 
'human'  vision  can  recognize  partial  shapes  with  no 
extra  effort,  recognition  of  partial  shapes  should  be  a 
built-in  part  of  shape  analysis  algorithms.  The  specific 
reason  why  an  algorithm  using  only  curvature  as  a 
criterion  does  not  work  well  for  comparing  partial 


39 


shapes  can  be  explanied  as  follows: 

A  partial  shape  contains  a  borrowed  segment  which 
belongs  either  to  another  occluding  object,  or  to  the 
boundary  imposed  at  the  limits  of  the  field  of  view,  or 
to  itself  but  not  to  the  complete  shape.  In  other  words 
a  partial  shape  contains  one  or  more  segments  which 
belong  to  the  complete  shape  and  some  which  do  not.  The 
task  of  comparing  a  partial  shape  to  a  complete  one  is 
that  of  locating  segments  on  both  complete  and  partial 
shapes  respectively,  such  that  their  log  (  <  )  plots 
differ  by  an  additive  constant.  This  constant  can  be 
eliminated  by  subtracting  the  means  over  the  respective 
segments.  The  problem  with  this  solution  is  that  it  is 
always  possible  to  find  a  segment  of  length  * ,  on  the 
partial  shape,  (  where  i  can  be  made  as  small  as 
desired)  such  that  by  appropriate  interpolation  and 
resampling  within  the  limits  of  resolution,  the 
curvature  of  this  segment  can  be  made  to  match  the 
curvature  of  some  segment  on  complete  shape.  Therefore 
when  no  syntactic  knowledge  is  given,  or  the 
relationship  between  the  size  and  the  number  of  samples 
is  unknown,  it  is  not  possible  to  decide  on  the 
necessary  length  of  the  segments  to  be  matched,  before 
it  could  be  concluded  that  the  partial  shape  is  a  part 
of  the  complete  shape.  Bence  without  a  priori  knowledge 


about  the  complete  shape,  matching  a  few  segments  of  the 
partial  shape  to  a  complete  shape  would  yield 
inconclusive  results. 

It  is  imperative  therefore  to  derive  some  syntactic 
knowledge  about  the  shape.  For  some  shapes  such  as  the 
ones  shown  in  Fig.  2.,  where  the  peaks  and  valleys  in 
the  curvature  correspond  to  the  vertices  or  critical 
point  of  the  shape,  enought  syntactic  knowledge  can  be 
derived  by  finding  the  peaks  and  valleys  in  the 
curvature  plot  to  permit  a  comparison. 

But  for  shapes  like  the  one  shown  in  Fig.  3, 
sufficient  information  cannot  be  derived  from  the 
curvature  plot  alone  to  permit  a  comparision.  Consider 
the  curvature  plot  of  the  shape  shown  in  Fig.  5.  which 
is  shown  in  Fig.  9 .  The  plot  of  Fig.  9  .  which  was 
obtained  by  using  (III. 8)  ,  has  only  one  peak 
corresponding  to  the  apex  of  the  cardioid.  Also  note 
that  the  curvature  function  obtained  by  using  (m. 8) 
for  the  cardioid  before  and  after  rotation  are  the 
identical  as  predicted,  in  sharp  contrast  to  the 
differences  demonstrated  in  Fig.  6-a  and  6-b.  Thus  if 
only  the  peaks  and  valleys  in  the  curvature  were  used  to 
derive  the  syntactic  knowledge  about  the  shape,  then  the 
cardioid  would  be  described  by  only  one  point,  namely 


the  apex.  This  is  clearly  insufficient.  Therefore,  peaks 
and  valleys  derived  from  curvature  alone  cannot  be 
uniformly  used  as  a  single  entity  to  represent  a  shape. 
This  is  not  to  imply  that  curvature  is  not  a  useful 
parameter,  but  rather  to  demonstrate  that  additional 
information,  or  attributes,  are  needed  to  fully  describe 
the  shape. 


CHAPTER  IV 


IV. a  LINE  OP  SIGHT  CONCEPTS. 


The  set  of  points  which  define  a  shape  may  be 
considered  as  a  Fuzzy  set  in  which  various  points  are 
assigned  to  it  with  various  degrees  of  membership. 
Defining  a  precise  criteria  on  which  a  degree  of 
membership  can  be  assigned  to  a  point  on  a  shape  is  very 
difficult.  However  past  reseachers  have  discovered  [14] , 
[19]/  [32]/  [911,  that  the  points  which  should  have  a 
higher  degree  of  membership  than  the  rest  are:  1) 
Points  of  Maxima  2)  Points  of  Minima  3)  Points  of 
Inflection  4)  Points  of  Intersection  5)  End  points  of 
open  curves  6)  Points  where  the  curvature  changes  sign 
or  magnitude. 

In  the  past  shape  analysis  efforts,  points  of  maxima, 
minima,  and  inflection  points  of  curves  have  been 
extracted  from  the  shape  data  without  due  consideration 
of  the  coordinate  axes.  This  approach  inevitably  leads 
to  errors,  because  these  points  have  no  meaning  unless 
the  coordinate  axes  are  first  defined.  The  problem  with 
such  an  approach  is  that  if  a  set  of  coordinate  axes  is 
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chosen  independent  of  the  shape  under  analysis  then  the 
maxima,  minima,  and  points  of  inflection  will  not  be 
independent  of  rotation,  of  the  shape,  with  respect  to 
the  chosen  coordinate  axes.  The  conclusion  is  that  the 
coordinate  axes  upon  which  the  maxima,  minima,  and 
points  of  inflection  are  based  must  be  dependent  on  the 
shape  itself. 

The  next  logical  question  is  whether  there  should  be 
one  set  of  coordinate  axes  or  many?  The  answer  is  not 
straight  forward.  Some  shapes  require  more  than  one  set 
of  coordinate  axes  while  others  require  only  one.  Before 
discussing  the  method  for  determining  critical  points, 
two  definitions  are  presented  which  will  be  used  in  the 
sequel . 

Definition  I:  A  curve  c  is  said  to  be  in  Line  of  Sight 
of  a  point  P  (LSP)  if  every  point  on  the  curve  c  can  be 
connected  to  P  without  intersecting  the  curve  at  any 
other  point.  Otherwise  the  curve  is  said  to  be  Not  in 
Line  of  Sight  of  the  point  P  (NLSP) . 

Examples  of  (LSP)  curves  of  a  point  P,  where  P  is  the 
centroid  C  of  the  shape,  which  are  then  denoted  as 
(LSC) -curves  ,  are  shown  in  Fig.  10.  The  concept  of 
Line  of  Sight  of  a  Point  is  useful  not  only  in  analysis 
of  shapes  but,  it  also  helps  in  reducing  the  number  of 


computations  involved  in  the  ALS  method. 

Definition  II:  Line  of  Sight  of  a  Straight  line  axis  L: 
Let  n  be  the  normal  projection  of  a  curve  c  onto  a 
straight  line  L.  The  curve  c  is  then  said  to  be  in  Line 
of  Sight  of  Axis  L  (LSA)  if  all  points  from  c  can  be 
mapped  injectively  onto  n. 

Examples  of  curves  which  are  in  Line  of  Sight  of  a 
single  straight  line  Axis  (LSA)  are  shown  in  Fig.  11, 
along  with  curves  Not  in  Line  of  Sight  of  a  single 
straight  Axis (NLSA) . 

In  the  Adaptive  Line  of  Sight  method,  which  is 
discussed  next ,  the  shape  is  divided  into  a  set 
of  segments  which  are  in  line  of  sight  of  a  set  of  axes. 
Dividing  the  shape  into  this  set  of  segments  is 
equivalent  to  defining  the  shape  in  terms  of  single 
valued  functions.  It  follows  from  the  definition  of 
single  valued  function  that  fewer  ambiguities  should 
result,  that  is,  the  maximas,  minimas,  and  inflection 
points  are  now  determined  from  single  valued  functions 
rather  than  a  multivalued  function.  The  actual  method 
for  deterimining  critical  points  is  presented  in  the 
next  section. 


EXAMPLES  OP  SlNCLY  CLOSED  CURVES 


IV.  b  ADAPTIVE  LINE  OF  SIGHT  METHOD. 


In  the  Adaptive  Line  of  Sight  method  the  critical 
point  determination  is  based  on  a  set  of  coordinate  axes 
that  are  dependent  on  the  shape  under  consideration.  As 
previously  discussed,  this  will  allow  critical  points  to 
be  determined  with  fewer  ambiguities.  The  procedure  is 
adaptive  in  the  sense  that  it  adapts  to  the  shape  data 
under  consideration.  This  will  become  evident  in  the 
sequel . 

The  Adaptive  Line  of  Sight  method  is  a  two  pass 
process.  In  the  first  pass,  the  shape  is  divided  into  a 
minimum  number  of  segments,  or  parts,  by  an  appropriate 
set  of  'critical  points'.  The  members  of  this 
segmenting  set  of  critical  points  are  defined  to  be 
those  points,  such  that  all  boundary  points  in  between 
any  two  adjacent  critical  points  have  the  following  two 
properties  with  respect  to  the  straight  line  L  joining 
every  two  adjacent  critical  points: 

1)  the  boundary  points  are  on  the  same  side  of  straight 
line  L  joining  the  adjacent  critical  points, 

2)  and  the  points  are  in  line  of  sight  of  L. 


Prom  the  definition  of  line  of  sight  (Definition  II) 
it  is  clear  that  this  means  that  the  boundary  curve  has 
a  unique  one  to  one  projection  on  L.  It  should  be 
emphasized  that  the  minimum  number  of  critical  points 
are  obtained  during  the  first  pass.  This  is  done  by  an 
exhaustive  search  process  that  locates  all  points  such 
that  the  above  properties  are  satisfied.  Often  the 
minimum  is  not  unique ,  in  which  case,  the  required  set 
is  obtained  by  summing  all  the  minimal  sets. 

In  the  second  pass,  the  points  of  maxima,  minima  and 
inflection  between  every  two  critical  points  are 
detected  using  the  derivative  of  the  normal  distance  of 
the  point  from  L  with  respect  to  the  distance  along  L.  A 
moving  average  of  these  normal  distances  may  be  used  to 
eliminate  the  effect  of  noise.  The  critical  points  found 
in  pass  one,  the  segmenting  set,  and  the  critical  points 
found  in  pass  two,  are  defined  to  be  the  members  of  the 
Fuzzy  shape  set  with  the  the  degree  of  membership 
depending  on  the  number  of  minimal  sets  they  are  found 
in.  A  through  description  of  both  the  computational  and 
detection  aspects  of  the  Adaptive  Line  of  Sight 
Algorithm  along  with  a  complete  flowchart,  is  given  in 
the  Appendix. 

The  algorithm  was  used  on  the  shapes  shown  in  figures 


2-a,  2-b,  3,  12,  and  13.  The  shapes  in  figures  2-a  and 
2-b  were  generated  by  sampling  the  curve,  obtained  by 
linearly  interpolating  between  the  vertex  points,  at 
constant  intervals  of  arc  length.  The  number  of  samples 
for  these  shapes  was  arbitarily  chosen  to  be  138  and  273 
respectively.  The  cardioid  shown  in  Fig.  3.  which  was 
generated  by  sampling  the  function:  r  *  60  (  1  +  cose  ) 
at  constant  intervals  of  6  and  not  the  arc  length,  has 
256  samples.  The  shapes  shown  in  figures  7  and  8  were 
obtained  by  sampling  two  hand  sketched  shapes,  through  a 
tablet,  which  was  used  as  an  input  device,  for  the 
Tektronics  4081  graphics  system.  The  sampling  process 
being  physical,  resulted  in  non-uniform  sampling.  These 
shapes  have  200  and  280  samples  respectively.  No 
interpolation  or  resampling  was  carried  out  before  the 
data  representing  these  shapes  was  input  to  the  ALS 
algorithm. 


The  critical  points  obtained  using  the  Adaptive  Line 
of  Sight  of  algorithm  on  these  shapes  is  shown  in 
figures  14,  15,  and  16.  Some  of  these  shapes  were 
rotated  to  confirm  that  the  algorithm  was  independent  of 
rotation.  It  can  be  seen  from  these  figures  that  in  the 
absence  of  predefined  resolution  the  critical  points 
tend  to  occur  in  clusters,  necessitating  a  post 
processing  step,  such  as  replacing  the  clusters  with  a 
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PIG.  16'  .  THE  CRITICAL  POINTS  OP  TBS  NONONIPORMLT  SAMPLES  SHAPES 


TYPICAL  MINIMAL  SETS. 
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single  point 


Note  from  the  shape  shown  that  the  ccitical  points 
are  very  close  to  those  perceived  by  human  vision, 
despite  of  the  total  dissimilarity  in  the  shapes.  Also 
note  that  the  post  processing  is  done  twice.  First  time 
clusters  are  temporarily  replace  by  single  points  to 
count  the  number  of  critical  points.  The  final 
replacement  of  the  clusters  by  single  points  occurs 
when  all  the  minimal  sets  have  been  superimposed. 

In  the  next  chapter  the  basic  concepts  of  size 
variable,  shape  vector  and  shape  space  are  introduced 
they  are  then  used  in  chapter  V  for  organizing  critical 
points  into  feature  vectors  or  subshape  vectors. 


CHAPTER  V 


BASIC  SHAPE  CONCEPTS 


V.a  MEASUREMENT  VECTOR. 


In  this  chapter  several  general  concepts  from 
allometric  disciplines  [52],  are  combined  with  some 
entirely  new  concepts  concerning  shapes.  This 
combination  provides  the  foundation  for  a  new  approach 
to  defining  shape  vectors  in  the  appropriate  shape 
space.  This  new  space  is  a  combination  of  properties  of 
vector  spaces  and  is  called  definitively  shape  space. 

Applying  the  shape  space  concepts  to  the  shape 
analysis  problem  provides  a  basis  for  the  recognition 
that  the  features  of  two  or  more  shapes  under  analysis 
are  the  same.  For  instance,  with  shape  space  concepts  it 
is  possible  to  determine  that  a  partial  shape, 
independent  of  size  or  rotation,  belongs  to  a  more 
complex  whole  shape. 

Typically  in  a  shape  analysis  problem,  an  algorithm 
operates  on  the  shape  data  according  to  a  set  of 


criteria  for  the  purpose  of  reaching  a  decision  of  some 
sort.  Usually  the  decision  is  whether  or  not  two  or  more 
shapes  are  the  same.  The  shape  analysis  algorithm 
utilizes  measurements  such  as  curvature ,  and 
measurements  between  predefined  points  on  the  shape, 
such  as  length  width,  diameter,  area,  etc.  Therefore,  if 
K  shapes  are  under  analysis,  and  N  measurements  ma)|  are 
made  on  each  shape,  then  the  result  is  the  K  measurement 
vectors, 
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where  the  first  subscript  n  refers  to  the  nth 
measurement  between  the  predefined  points  on  the  k  th 
shape.  The  shape  or  the  object  is  represented  by  the 
second  subscript.  All  measurements  are  assumed  to  be 
positive  and  the  angle  is  the  angle  made  by  the 

nth  measurement  on  the  kth  shape  with  respect  to  a 
reference  direction.  Any  two  objects  will  then  be  said 
to  have  the  same  shape  with  respect  to  these 
measurements  if  one  vector  is  a  scalar  multiple  of  the 
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Mj  -  a  (V. 2) 

where  a  is  a  scalar  greater  than  zero. 


V.b  SIZE  VARIABLE 

The  geometric  significance  of  (V.2)  is  that  in  the  N 
dimensional  space  of  positive  measurements  all  points  on 
a  straight  line  through  the  origin  define  the  same 
shape.  Points  can  be  uniquely  located  on  the  positively 
directed  line  by  finding  the  intersection  of  first  order 
surfaces  with  the  line  defining  the  shape.  The  class  of 
functions  which  define  these  surfaces  are  homogenous 
functions  of  order  one,  of  the  measurements,  m|k  n»l, 
2,.N.  The  mathematical  representating  for  this  class 
is, 

*  (  »*)  *  *  5  (  ■»„) 

v  (V. 3) 

2  (  )  »  0 

where  a  is  scalar 

and  S  refers  to  a  countably  infinite  class  of 


functions,  with  members  Z-,  (  ) 


The  distance  from  the  origin  to  the  intersection  of  a 
particular  member  of  the  class  2j  (  m, k  )  with  the  ray 
through  origin  defining  the  shapef is  refered  to  as  the 
size,  scale  factor  or  the  normalization  factor  with 
respect  to  that  member.  Following  Mosaiman  [52] 
terminilogy,  in  the  sequel  Zj  (  )  will  be  refered 
to  as  the  size  variable.  Some  examples  of  size  variable 
in  a  measurement  space  of  two  vectors  are  shown  in  Fig. 
17a  and  b. 

When  comparing  shapes,  it  is  necessary  to  choose  a 
size  variable  that  is  independent  of  the  shape  .  under 
comparison.  A  necessary  requirement  for  the  verification 
a?,  this  statement  is  the  defintion  of  a  shape  vector 
[52] ,  as  a  vector  valued  function  of  vectors,  and  the 
concept  of  shape  space. 

In  all  that  follows  it  will  be  assumed  that  the 
measurements  are  always  made  between  the  centroid  of  the 
shape  in  question  and  some  other  fixed  point  on  the 
shape. 
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PIG.  17  .  EXAMPLES  OP  SIZE-VARIABLE  Dl  2-0  SHAPE  SPACE 


V.c  SHAPE  VECTOR 


The  Shape  vector  Skl  is  defined  to  be  the  ratio  of  the 
measurement  vector  Mk  to  the  size  variable  Zj. 

* "  ■  "»/  zi 

S^^f  m^^lk  ,  .  .  .  ) 

Zj  (  )  Z  j{  )  Zj(  )  (V.  4) 

It  may  be  noted  that  the  first  superscript  of  S*1 
corresponds  to  the  object  whose  shape  is  under 
consideration  while  the  second  superscript  corresponds  to 
the  size  variable  chosen  from  the  class. 

Here  the  points  to  be  emphasized  are  that, 

1 )  All  measurements  are  made  between  pre-defined  points. 

2)  A  shape  vector  is  defined  with  respect  to  a  size 
variable  Z  j  (  m^ )  .  Thus  only  shape  vectors  defined 
with  respect  to  the  same  size  variable  can  be  compared. 

3)  If  two  shape  vectors  are  equal,  then  the  two  objects 
have  the  same  shape  with  respect  to  the  measurements. 

4)  The  shape  vector  should  be  independent  of  the  size 
variable. 


These  statements  merit  further  discussion  and 

clarification  because  of  their  implications.  Consider 
for  example  the  two  different  measurement  vectors 
extracted  from  the  simple  shape  shown  in  Fig.  18.  The 

shape  is  a  unit  square.  In  the  first  case,  the 

measurement  vector  is, 

nijj  /  »21  *  '  ®4i  )*(  1  ,1 .414,  1) 

(V.5) 

while  in  the*  second  case,  the  measurement  vector  is 
Mj'-C  m^  ,m22  '  m32  m42  >-(  lr  lf  lf  1}  • 

(V.6) 

where  the  '  in  the  above  equation  indicates  that  the 
measurements  are  between  a  different  set  of  feature 
points  of  the  shape.  Observe  that  in  the  above  two 
equations,  the  measurement  vectors  are  considered  as  a 
function  of  the  distances  only.  Comparison  of 
measurement  vector  as  a  function  of  the  angle  a#|,  is 
more  complex  and  is  delayed  till  chapter  V.  Now  if, 

Zj(  mik  )  -  m3k  (V.7) 

is  chosen  as  the  size  variable  then  the  corresponding 
shape  vectors  are. 


S-C.707  ,  .707  ,  1  ,  .707) 


(V.  8) 


S2-(  1  ,  1,1,  1)  (V. 9) 

Comparing  these  two  shapes  without  any  reference  to  the 
size  variable  or  the  points  between  which  the 
measurements  were  made/  one  would  conclude  that  the  the 
two  shapes  are  not  the  same.  Thus  not  only  has  the 
functional  form  of  the  size  variable  to  be  the  same  but 
the  measurements  involved  in  the  functional  relationship 
have  to  be  between  the  same  feature  points.  Now  consider 
the  shapes  in  Fig.  19.  This  is  the  typical  situation  in 
which  occlusion  occurs  or  the  shape  is  outside  the  field 
of  view  of  a  camera  (or  some  other  measuring  device) . 
Assume  that  both  shape  boundaries  are  represented  by  an 
equal  number  of  samples  and  that  every  point  on  each 
shape  boundary  is  defined  as  a  feature  point.  Now,  if  the 
standard  deviation  of  the  data  is  chosen  as  the  size 
variable  then  it  is  obvious  that  the  standard  deviation 
of  the  two  shape  boundaries  are  different.  Therefore/ 
this  is  a  case  where  the  size  variable  is  dependent  on 
the  shape/  which  implies  that  the  shape  vector  is 
dependent  on  the. size  variable.  Comparision  of  the  two 
shape  vectors  with  such  a  size  variable  is  bound  to  lead 
to  errors.  These  problems  can  be  alleviated  by  defining 
the  quantities  in  the  proper  context.  This  can  be 
accomplished  only  if  the  variables  are  defined  in  the 


proper  space. 


V.d  SHAPE  SPACE 


The  space  is  defined  in  terms  of  its  properties  in 

the  usual  manner  and  then  two  theorems  addressing  the 

problem  of  partial  shapes  are  stated  and  proved. 

Assuming  that  Sk*  {  sj1  )  is  a  shape  vector  consisting 

of  well  defined  operations  on  the  measurements  of  the 

shape  under  consideration.  These  measurements  are  the  m| 

ki  ki 

' s  previously  defined;  the  elements  of  s  (  )  are 

obtained  by  the  following  operation, 

si'  -  alh^.k  /Z.(  mik  )  (V.10) 

The  shape  vector  Sh*  (  S*J)  must  satisfy  the  following 
properties  in  addition  to  the  properties  of  normal 
Euclidean  space. 

PROPERTIES  OF  SHAPE  SPACE 

1)  The  shape  vector  is  independent  of  the  size  variable. 
This  implies  that, 

Skl  (  a  sj*  )  •  a  Skl  (  sji  )  (V.ll) 

where  n  »  1,  2, . ,N,  k»l,  2,  ....  K, 

and  a  is  a  scalar. 

2)  The  shape  vector  is  independent  of  translation  that  is 

sli  (•{*♦  •,  )  -  »*'  (."  ) 

where  n  •  1,  2  ,  . N,  K»l,  2, 


(V.  12 


(V. 13) 


and  Sgis  a  constant  vector. 

3) The  shape  is  independent  of  rotation. 

*"  (  U,  1  .  a“  1  .(|1  1 
where  n«l,  2,  . . . . .N  k»l,  2,  ....  K, 
and  U|  is  a  constant  angle. 

The  vector  obtained  by  using  a  set  of  measurements  on 
a  partial  shape  must  still  be  contained  in  the 
space. Unless  a  size  variable  is  found  which  is 
independent  of  both  shapes  (  both  the  partial  and  the 
complete)  it  is  not  meaningf ull  to  compare  the  shape 
vector  in  shape  space.  It  is  not  possible  to  find  a 
size  variable  which  is  a  totally  independent  continious 
function  of  measurements  made  on  both  shapes.  This  being 
the  case  the  only  choice  left  is  to  split  the  shape  into 
parts  or  subshapes  and  define  a  size  variable  which  is 
piecewise  continuous  over  these  parts.  Two  theorems 
relating  the  subshape  to  the  complete  shape  in  shape 
space  are  stated  and  proved.  These  theorems  are  used 
extensively  in  the  sequel. 


V.e  THEOREMS  ON  CONCATENATED  SHAPE  VECTORS. 


THEOREM  1  t  The  vector  formed  by  concatenating  a  series 
of  shape  vector  is  a  shape  vector. 


Proof:  It  is  required  to  show  that  the  shape  vector 
resulting  from  concatenating  a  series  of  shape  vectors 
satisfy  the  three  properties  of  shape  space. Let  the 


concatenated  shape  vector  be 

3ei  -(S11  (  s“  ),  S22  (  ),  .  ., 


S  (  s  )  »  «  • ) 


where  u  ■  1  . . .  0,  v  •  1  . . .  V,  and  w  »  1  . . .  w.  and 

the  superscript  c  on  S  indicates  that  it  is  obtained  by 

ci  t  j 

concatenating  other  vectors.  Each  element  of  S  (  s, 
can  be  represented  by 


*■  *  m«e  /zi  <  ®«e  >  (V*1 
j  ®  1  ...  J?  and  n  ■  1  ...  N 

where  J  is  the  total  number  of  concatenated  vectors, 
while  N  is  equal  to  the  sum  of  the  total  number  of 
components  of  all  the  concatenated  vectors. 


Therefore, 


(8 1  )  “(  »  Sj  •  *3  ' 


,  sa  ,  . . .  • ) 


Multiplying  each  component  of  (V.14)  by  the  scalar  a 


(V.  1 8 


but  from  the  {  V.10)  we  have, 

Shfc  (  a  sj1  )  -  a  Shfc  (  sjh  ) 

using  (V.17)  in  (V.18)  results  in 
a  s^  )»(a  sj1  ,  a  s^  , .  .  . ,  a  sj*  ,  .  .)  (V.19) 

or 

a  S‘>  (  a,'<  )  .  sei(  a  aji  ). 

Equation  (V.20)  proves  property  i) 

To  prove  that  the  shape  vector  satisfies  proverty  ii)  it 
is  only  necessary  to  observe  that  each  member  of  the 
concatenated  vector  is  a  shape  vector .Therefore  each 
sf  ies: 

sCi  (  *«i  +  v sCi  (  ]  (V.21) 

The  proof  for  property  iii)  follows  in  a  similar  manner 
Q.E.D 

Theorem 2 :  The  shape  defined  by  a  shape  vector  obtained 
by  concatenating  a  series  of  shape  vectors  is  unique  if 
and  only  if  the  size  variable  and  size  defining  each 
member  of  the  set  are  known. 

Proof:  — .  The  surface  of  a  homogeneous  function  of 


order  one  will  intersect  a  positive  directed  straight 


line  at  only  one  point.  Since  a  shape  is  defined  as  a 
point  on  the  shape  ray,  if  follows  that  any  point  on 
this  ray  can  be  uniquely  determined  by  its  intersection 
of  a  size  variable  which  is  defined  as  a  homogenous 
function  over  the  positive  quadrant. 

:< - .  Assume  otherwise.  Then  there  is  at  least  one 

)  shape  vector  in  the  concatenated  set  whose  size 
variable  Zk(mwk)  can  be  choosen  arbitarily.  Then  the 
concatenated  vector  under  this  assumption  would  still 


satisfy  the  properties  of  the  shape  space,  i.e., 


a  SeJ  (  s^  ) 


S*1  (a  s',*) 


(V. 22) 


Now  the  R.H.S  of  (  V.22)  can  also  be  expanded  as 
-(as11  ,  a  S22  ,  ...  ,a  sfcfc  ....) 


(V. 23) 


(  _aMi. 

Zl 


•  f  •  '  ) 

zk  rv  ■ 


(V. 24) 


but  since  can  be  choosen  arbitrarily  as  long  as  it 

satisfies  the  definition  of  a  size  variable(  v.3). 


Choose 


zi  (  a  mB|  ) 


(V.25) 


substituting  (V.3)  and  (V.25)  if  (V.24)  the  following 
relationship  is  obtained 


CHAPTER  VI. 


FEATURE  SELECTION  &  COGNITIVE  STEP 


VI. a  FEATURE  VECTOR  FORMATION  PROCESS 


The  features  of  a  shape  are  essential  to  defining  the 
shape  in  terms  of  parameters  that  can  be  ultimately  used 
-  by  machine  for  decision  purposes.  However,  the  manner  in 
which  the  feature  defining  procedure  can  be  selected  is 
quite  variable.  Since  a  dependable  feature  selection 
procedure  is  fundamental  to  the  shape  recognition 
problem,  it  is  essential  that  this  aspect  of  shape 
recognition  be  addressed  with  specificity.  This  point  is 
punctuated  when  it  is  realized  that,  irrespective  of  the 
method  of  defining  and  detecting  critical  points,  a 
cognitive  algorithm  is  still  required  which  examines  in 
some  sense,  the  critical  points  of  the  shape  for  the 
purpose  of  reaching  a  decision  about  some  aspect  of  the 
shape.  Consider,  for  example,  the  shape  shown  in  Fig. 
2-a,  a  shape  such  as  this  swept  wing  plane  may  have 
thirty  to  fifty  critical  points.  The  human  eye  makes 

numerous  measurements,  automatically  and  sub- 
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consciously,  between  the  feature  points  and  determines 
their  relationship  to  one  another.  The  'most  important' 
— waofca,  with  the  relationship 

between  them,  comprise  the  decision  set.  The  term  ”  most 
important  "  is  difficult  to  define  mathematically, 
because  it  is  the  result  of  training.  An  unrefined 
cognitive  procedure  must  therefore  consider  the  set  of 
all  possible  measurements  between  the  critical  points. 
Obviously  this  is  a  very  large  number  of  measurements 
even  for  numbers  as  modest  as  thirty  to  fifty.  It  is 
well  known  that  this  totality  of  measurements  between 
critical  points  is  not  essential  to  the  decision 
process. 


It  is  necessary  therefore  to  determine  methods  for 
acquiring  the  minimal  set  of  measurements  or  features 
required  for  the  decision  process.  The  human,  apparently 
places  heavy  weighting  on  features  that  are  formed  by 
critical  points  that  are  symmetrically  opposite  about  an 
axis  and  features  that  are  extracted  from  adjacent 
critical  points  concerning  the  shape.  Without  any  prior 
knowledge  a  human  can  find  the  sets  of  axes  about  which 
some  critical  points  are  symmetrically  placed  with  very 
little  effort.  However,  such  a  task  is  almost 


insurmountable  for  a  machine  based  algorithm  unless  it 
is  performed  at  a  post-cognitive  level.  In  the  absence 


of  noise,  machine  recognition  (  cognitive  )  algorithms 
perform  reasonably  well  by  using  only  features 
^mg^tin^ofadjacentcritical  points. 

A  cognitive  algorithm  that  utilizes  measurements  such 
as  these  in  a  continuous  sequential  manner  would  be 
entirely  adequate  if  the  algorithm  for  detecting 
critical  points  is  totally  immune  to  noise,  round  off, 
and  truncation  errors.  For  example  any  extra  critical 
points  that  are  the  result  of  a  burst  of  noise  would 
prevent  any  continuous  sequential  recognition  algorithm 
from  yielding  conclusive  results. 

One  manner  by  which  this  problem  can  be  circumvented 
is  to  divide  each  of  the  shapes  under  analysis  into 
subshapes  in  terms  of  their  features,  and  then  compare 
the  features  of  these  subshapes  and  the  manner  in  which 
they  are  related  to  each  other.  It  is  recalled  from 
chapter  V  that  the  properties  of  shape  space  dictate 
that  the  measurements  which  define  a  feature  must  be 
made  from  the  centroid  of  the  set  of  critical 
points.  It  is  necessary  then  to  determine,  in  some 
manner,  the  minimum  number  of  critical  points  that 
adequately  define  a  feature.  If  each  feature  were 


defined  in  terms  of  only  two  critical  points,  then  all 
features  would  be  identical  since  the  comparisons  are 


made  with  respect  to  the  same  shape  independent  size 
variable.  Therefore,  the  minimum  number  of  critical 
points  that  can  form  a  distinguishable  feature  is  three, 
and  these  must  be  adjacent.  However  unless  the 
relationship  between  these  three  point  features  with  its 
adjacent  features  is  also,  considered  any  comparison 
(cognitive  )  algorithm  almost  always  leads  to 
ambiguities.  The  reason  for  this  is  because  the  three 
point  feature  forms  a  triangle.  An  examination  of  the 
shape  shown  in  Fig.  2-a  shows  that  it  contains  many 
similar  triangles. 

Unfortunately  the  mathematics  required  to  obtain  the 
optimal  number  of  critical  points  that  should  form  a 
feature  is  not  yet  developed.  Therefore,  it  is  necessary 
to  resort  to  the  psychological  aspects  of  the  human 
recognition  and  decision  process  as  well  as  the 
practical  aspects  such  as  the  implementation  and 
computational  requirements.  These  criteria  lead  to  the 
features  being  selected  as  follows: 

1)  Reconstruct  the  shape  by  connecting  all  the  adjacent 
points  by  a  straight  line.  This  is  called  the  critical 
shape  boundary. 

2)  The  feature  F^  is  then  formed  by  including  critical 
points,  Ct  and  at  most  three  adjacent  critical  points  on 
each  side  of  Cl  .  The  critical  points  chosen  must  be  in 


line  of  sight  of  C|  .  This  means  that  it  must  be 
possible  to  draw  a  straight  line  form  C|  to  each  of  the 
other  critical  points  defining  the  feature  without 
intersecting  the  critical  shape  boundary.  The  feature 
obtained  by  this  procedure  cannot  be  defended  as  optimal 
in  any  mathematical  sense.  However,  it  correlates  quite 
well  with  those  entities  that  humans  consider  features. 
The  features  corresponding  to  critical  points  and 

C  are  shown  in  Fig.  20. 
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A  desirable  improvement  to  the  above  feature  defining 
procedure  is  an  algorithm  for  deciding  whether  the  line 
joining  Ch  to  another  critical  point  in  the  feature  lies 
inside  or  outside  the  shape. 

The  cognitive  step  requires,  as  usual  a  dictionary  of 
the  features  of  the  complete  shape  against  which  the 
partial  shapes  are  to  be  compared.  The  partial  shape 
dictionary  will  henceforth  be  referred  to  as  the  problem 
text.  One  page  of  the  complete  shape  dictionary  is  shown 
in  Table  1.  This  page  contains  features  twenty-one  and 
twenty-two  of  the  swept  wing  plane  of  Fig.  21.  The 
table  includes,  in  addition  to  the  feature  number,  the 
critical  points  of  that  feature  along  with  their  x  and  y 
location,  the  x  and  y  location  of  the  centroid  of  all 
the  critical  points  contained  in  the  feature,  the  size 


FIG.  21  .  CRITICAL  POINTS  AND  EXAMPLES  OF  FEATURES  FOR  THE  PLANE. 


of  the  feature,  the  normalized  components  of  the  shape 
vector,  the  angle  of  the  shape  vector  component  and  the 


of sight  code.  The  of sight  code  is  a  binary  code  of 
length  seven,  associated  with  with  the  kth  feature  k  «  1 
.  .  .K.  The  nth  bit  of  the  code  is  equal  to  one  if  the 
critical  point  Cj_b  is  in  line  of  sight  of  the 
critical  point  Ck,  where  I  »k+4  and  n  -1...7. 


The  normalized  shape  vector  component  are  defined  as 


s.ki- 


x.2 


-.2 


sqrt  ( (x,h  -xj  +  (y>k  -  y)  )/  Zj 


(VI. 1) 


where  jthe  size  variable  was  chosen  to  have  the  form, 

{  (x„k-  x)2  +  (yfk-  y)2  }  (VI. 2) 

where  'the  subscript  k  refers  to  the  feature  number. 


Zj  - 


while  the  subscript  n  refers  to  the  measurement.  The 
angle  alpha  is  measured  with  respect  to  a  fixed 


reference,  and  is  obtained  by  the  following  equation, 
'C 


L 


“ik  * 


arctan{(y<k  -  y)/(x„k  -  x) } 
if  (aft  -  x)  >  0 
(*/2)  +  arctan{(y|k  -  yl/fx,!,  -  x)} 
if  (x,(  -  x)  <  0 


(t/4)  if  ( X|k  -  x)  ■  0  and  (y>k  -  y)  > 


_C3» /4)  if  (X,),-  -  x)  ■  0  and  (y#k  -  y)  <  0 

(VI. 3) 

These  same  quantities  are  obtained  to  form  a  dictionary 
for  the  partial  shape  (problem  text) .  Each  page  of  both 


the  dictionaries  begins  with  a  feature  set  that  is  not  a 
subset  of  another  feature  set.  This  feature  set  is 
defined  as  the  uncovered  feature  set.  The  covered 
feature  sets  are  arranged  in  the  order  of  cardinality 
below  the  uncovered  feature  set  on  each  page.  The 
purpose  of  this  architecture  is  to  simplyfy  the 
computational  requirements  for  the  cognitive  step. 

In  general  all  the  features  in  the  dictionary  will 
not  be  contained  in  the  problem  text.  It  is  also  true 
that  the  problem  text  contains  features  that  are  not 
present  in  the  dictionary.  This  becomes  apparent  by 
examining  Fig.  20.  Therefore,  the  fact  that  a  feature 
is  contained  in  the  problem  text,  does  not  imply  that 
the  partial  shape  is  not  a  part  of  the  complete  shape, 
because  it  is  not  necessary  for  the  partial  shape  to 
have  fewer  points  that  the  whole  shape.  Therefore, 
further  examination  is  required  before  a  decision 
regarding  the  problem  text  can  be  made. 


VI. b  THE  COGNITIVE  PROCESS 

The  decision  procedure  consists  of  selecting  an 
arbitrary  word  from  the  problem  text  dictionary.  A 


problem  text  word  is  o£  coarse  a  feature  from  the 


partial  shape  under  comparison.  The  shape  vector  from 
the  problem  text  is  compared  to  the  shape  vector  in  the 
dictionary  always  starting  on  page  one  of  the 
dictionary.  The  comparison  continues  until  a  match 


occurs. 


VI. b. i  THE  METRIC  DISTANCE 


The  process  of  matching  distances  is  straight  forward 
and  can  be  achieved  using  linear  correlation  [23] , 
metric  difference,  or  any  other  metric.  The  following 
metric  was  used  here. 


SOMABSDIF 


(VI. 4) 


'  ’<  di 

ABS  {(S,  -  sj*  )} 


Vl.b.ii  THE  ANGULAR  CORRELATION 


The  angle  is  an  unreliable  variable  for  direct 
comparision.  Since  the  angle  is  circular  variable  [4] , 
methods  of  circular  statistics  should  be  used  for  it 
comparison.  To  illustrate  the  issue  at  hand  consider  a 


feature  vector  whose  angle  vector 
following  three  components. 


If  '  0.785 


has  the 


2P  ■  1.570  radians 


5 


Now^  if  the  feature  undergoes  an  angular  rotation  of  2 


radians  then  the  corresponding  angle  vector  aR^  will 
have  the  following  components  in  the  new  position. 


mm  m 

—  * 

“14 

2. 785 

“24 

m 

3.570 

“3*4 

0.720 

_  -4 

_  — 

radl ans 


The  last  component  is  equal  to  .72  radians  as  all 
operations  are  performed  modulo  6.28.  Mote  a  direct 

comparison  of  the  above  two  vectors  would  not  be 
conclusive.  Mow,  consider  the  angular  displacemet  vector 
between  »  and  °a4«  It's  components  are, 


radians. 


ft*1 

r 

2.00 

m 

2.00 

kJ 

-4.28 

Observe  that  even  though  a  constant  was  added  to  all  the 
component  of  aRj  to  obtain  a<a  the  difference  between 
“«»  and  aR-  is  note  a  constant.  Hence  one  might  again  be 
lead  to  believe  that  the  two  angle  vectors  are  not  the 
same.  Using  the  methods  of  circular  statistics  [4],  such 
pitfalls  can  be  easily  avoided.  One  way  of  testing 
whether  the  4^  ' s  from  which  the  sample  (components)  are 
drawn  are  unidirectional  or  whether  there  is  any 
statistical  evidence  of  directedness  is  by  the  using  the 
mean  vector  length  rasa  measure  of  correlation. The  mean 
vector  is  given  by, 


where  et*s  are  unit  vectors  each  pointing  in  the 
direction  of  the  1 s.  Now  if  the  resulting  vector 
length  is  denoted  by  R  then, 


n 

z*. 


Each  of  the  eB's  can  be  expressed  as, 

U  N 

e„  -  cos  4,  +  sin «„ 


(VI. 6) 


(VI. 7) 


(VI. 8) 


substituting  (VI. 8)  into  (VI. 9)  the  following  expression  for 


r  can  be  obtained 


»  1  {  (  ^  cos  i[*  )2  +  (  ^  sin  6,'  )2  } 


M  2  1/2 


(VI. 9) 


The  following  properties  sake  r  an  useful  parameter 
for  comparision  of  directedness  of  two  circular  vectors. 

1)  The  value  of  r  is  independent  of  the  zero 
■direction. 

2)  r  can  vary  between  0  and  1.  If  r  is  equal  to  1  then 
all  the  circular  variables  (  in  this  case  the  angular 
displacement  vector  components)  point  in  the  same 
direction,  while  r  is  equal  to  zero  indicates  that  the 

components  are  randomly  distributed.  Values  in  between  0 
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and  1  indicate  various  degrees  of  randomness  in  the 
directedness  of  the  components.  In  general  for  a  word 
to  match  a  feature  the  following  three  conditions  must 
hold  simultaneously, 

1)  The  SUMABSDIF  should  be  below  a  threshold 

2)  The  RVALUE  (  value  of  r)  should  be  above  a  threshold 

3)  The  word  and  the  feature  should  have  the  same  ofsight 
code. 


Vl.b.iii  LOCATION  OF  MISMATCHED  POINTS 

The  next  step  is  to  compare  the  next  problem  text 
feature  vector  in  order  of  cardinality,  to  the  next 
feature  in  the  dictionary  and  so  forth.  An  example  of 
this  technique  is  given  by  comparing  Tables  1  and  2.  In 
this  experiment  a  feature  vector  from  the  partial  shape 
was  selected  for  comparison.  It  should  be  emphasized 
that  the  feature  vector  is  from  the  problem  text  of  the 
partial  swept  wing  aircraft  shown  in  Fig.  20.  The 
partial  shape  has  been  rotated  and  shifted  as  well  as 
scaled  to  insure  that  any  direct  template  matching 
procedure  will  fail.  This  also  demonstrates  that  the 
concatenated  feature  vector  matching  procedure  described 
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here  is  independent  of  rotation,  size  and  location. 
Feature  vector  twelve  (  word  }  was  arbitrarily  selected 
from  the  problem  text.  By  comparing  tables  1  and  2  From 
of  the  plot  the  RVALUE  and  the  SUKABSDIF  shown  in  Fig. 
22-a  and  Fig.  22-b  it  is  apparent  that  word  12  matched 
feature  vector  21  of  the  dictionary.  It  should  be  noted 
that  this  match  occurs  even  though  the  location  of  the 
critical  points  and  the  centroid  of  the  feature  of  the 
partial  shape  are  different  from  those  same  quantities 
for  the  whole  shape  because  of  the  rotation  and  shift. 
In  this  way  a  correspondence  table  is  then  established 
between  the  critical  points  of  the  features  in  the 
dictionary  to  the  critical  points  of  the  word  in  the 
problem  text.  The  next  step  is  to  proceed  in  sequential 
order  to  the  next  word  in  the  problem  text  e.g.  word  13 
(n+1)  which  is  sequentially  next  to  word.  12  (n)  of  the 
problem  text  does  not  match  feature  vector  22  (m+1)  of 
the  dictionary.  However  since  feature  vector  22  is  on 
the  same  page  as  feature  vector  21  (  m  )  of  the 
dictionary  ,  word  13  is  mismatched  to  feature  22  because 
it  contains  a  critical  point  which  is  not  contained  in 
word  12.  An  examination  of  tables  1  and  2  indicates 
that  the  critical  point  Cj  is  contained  in  word  13  but 
not  in  word  12. 


The  mismatched  critical  point  is  first  compared  to 
the  correspondence  table.  Zf  it  is  not  found  in  the 
correspondence  table  then  it  is  stored  in  a  mismatch 
table.  At  any  latter  stage  a  mismatched  critical  point 
is  erased  from  the  mismatched  table  if  some  word 
containing  the  mismatched  critical  point  matches  some 
feature  of  the  dictionary. 

In  general  if  word  n  matches  feature  m  then  it  is 
expected  that  word  (  n  +1  )  will  match  feature  (  m  +  1) . 
If  feature  (m+1)  is  on  the  same  page  as  feature  m  then 
it  is  easy  to  isolate  the  mismatched  point  as  in  the 
above  example.  If  feature  a+1  is  not  on  the  same  page  as 
feature  m  or  the  concept  of  pages  is  not  used  then  in 
order  to  isolate  the  mismatched  critical  point  then  the 
feature  m+1  and  the  word  n+1  have  to  be  revised  into 
concatened  shape  vectors  with  each  subshape  vector  of 
three  measurements.  A  revised  feature  vector  for 
feature  22  is  shown  in  Table  3.  While  the  revised  words 
for  word  13  are  shown  Table  4.  A  comparison  of  the  above 
two  tables  again  isolates  Cj  as  the  mismatched  critical 
point.  It  must  be  noted  again  from  these  tables  that 
angle  being  a  multivalued  function  cannot  be  relied  upon 
as  a  variable  that  can  be  used  in  the  direct  comparison 
process. 
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The  usefulness  of  the  Fuzzy  shape  concept  is 
demonstrated  by  the  results  of  the  ALS  algorithm  on  the 
pig  of  size  B  and  the  partial  pig  of  size  10xB.  More 
than  one  minimal  set  can  be  obtained  by  starting  at 
various  points  on  the  back  of  the  pig  or  the  earlobe  of 
the  pig.  However, the  points  which  are  found  most  often 
in  minimal  sets  have  been  enclosed  by  circles  and  are 
shown  in  Fig  23-a  and  23-b.  These  clusters  or  points  of 
the  highest  degree  of  membership  are  retained  by  the 
algorithm  while  the  others  are  discarded.  The  effect  of 
postprocessing  on  the  clusters  of  highest  degree  of 
membership  gives  the  final  minimal  set.  The  final 
minimal  set  for  the  pig  and  the  partial  pig  are  shown  in 
Figures  24-a  and  24-b  respectively.  Word  21  of  the 
partial  pig  is  listed  in  table  6.  The  plot  of  the 
RVALUE  and  the  SUMABSDIF  of  word  21  when  compared  with 
features  with  the  same  ofsight  code  in  the  dictionary  is 
given  if  Fig  25-a  and  25-b.  From  these  plots  it  is  clear 
that  word  21  matched  feature  31.  The  details  of  feature 
31  are  listed  in  table  5. 
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CHAPTER  VII. 


CONCLUSION  AND  DISCUSSION 

VII. a  CONCLUSION 


A  system  for  the  recognition  of  shapes  which  is  in 
some  sense  similar  to  human  vision  system  is  described. 
The  system  is  called  the  partial  shape  recognition 
system  since  it  is  capable  of  recognizing  shapes  based 
on  comparison  of  properties  of  parts  of  a  shape  rather 
than  only  the  global  properties.  Though  not  all  aspects 
involved  in  implementing  the  system  were  discussed,  some 
of  the  necessary  aspects  were  detailed. 

A  new  concept  of  treating  shapes  as  vectors  in  shape 
space  is  introduced  and  described.  Also  two  theorems 
relating  to  the  process  of  comparing  partial  shapes  to 
the  complete  shape  were  stated  and  proved. 

A  new  procedure  of  determining  (  the  correspondence 
tokens)  the  critical  points  of  a  shape  is  described.  The 
procedure  is  named  the  Adaptive  Line  of  Sight  method.  In 
the  Adaptive  Line  of  Sight  method  the  critical  point 
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determination  is  based  on  a  set  of  coordinate  axes  that 
are  dependent  on  the  shape  being  examined.  Examples  were 
given  that  demonstrate  that  the  procedure  produces 
critical  points  that  are  independent  of  rotation,  size, 
displacement  and  correspond  closely  to  those  produced  by 
normal  human  cognitive  process. 

It  was  shown  that  the  critical  points  could  be 
organized  to  form  feature  vectors  using  the  Line  of 
Sight  of  a  Point  concept.  A  technique  for  comparing  the 
feature  vectors  of  a  set  of  shapes  is  described.  The 
comparison  procedure  is  based  on  syntactic  technique 
which  will  point  whether  one  shape  is  part  of  a  more 
complex  whole  shape,  or  whether  the  shapes  are  totally 
dissimilar. 

VII. b  DISCUSSION 

VII. b. i  CAN  COMPARISON  OF  SHAPES  BE  BASED  ON 
CRITICAL  POINTS  ALONE  ? 


Examples  presented  in  the  earlier  sections  showed 
that  the  shape  recognition  technique  based  on  comparing 


critical  points  of  the  highest  degree  of  membership  was 
sufficient  and  gave  good  results.  There  are  instances 
however,  when  a  natch  based  on  comparison  of  critical 
points  alone  is  insufficient.  A  simple  example  is  that 
of  an  ellipse.  The  critical  points  of  the  highest  degree 
of  membership  of  the  ellipse  shape  are  the  points  where 
the  ellipse  intersects  it  minor  and  major  axes.  These 
critical  points  are  also  the  critical  points  of  a 
polygon  (  which  looks  different  than  an  ellipse)  formed 
by  joining  the  adjacent  critical  points  by  a  straight 
line.  A  comparison  based  on  critical  points  alone  would 
not  be  able  to  make  out  the  difference  between  the 
polygon  and  the  ellipse.  In  general  the  shape 
represented  by  the  critical  shape  boundary  and  the  shape 
boundary  are  indistinguishable.  Another  example  is  that 
of  the  circle  and,  a  circular  arc  which  is  greater  than 
a  semi-circle  but  less  than  the  circle.  For  the  circle 
all  points  will  be  determined  to  be  equally  important. 
For  the  circular  arc  all  points  excluding  the  points 
shown  in  Fig.  26.  will  be  determined  to  have  the  same 
degree  of  membership.  The  points  shown  on  the  circular 
arc  have  a  higher  degree  of  membership  than  the  rest. 
Here  a  comparison  based  on  matching  critical  points  of 
the  highest  degree  of  membership  would  lead  to  wrong 
results.  In  such  instances  it  is  proposed  to  match  the 


’  AXIS  OF  SYMMETRY 


NOTE: ALL  POINTS  ON  THE  ARC  ARE  CRITICAL  POINTS  BY  THE  ADAPTIVE 
LINE  OF  SIGHT  METHOD;  HOWEVER  THE  POINTS  WITH  THE 
HIGHEST  DEGREE  OF  MEMBERSHIP  HAVE  BEEN  ENCLOSED  BY  CIRCULAR 
DOTS 

FIG.  26.  A  CASE  FOR  RELOCATION  OF  CRITICAL  POINTS  ? 


shapes  based  on  the  space  invariant  properties  o£  the 
segments  (  interpolated  and  resampled  if  necessary) 
between  critical  points.  Figure  1.  thus  shows  a  block 
"length  of  the  segments  to  be  matched”  which  is  also 
loaded  into  the  shape  data  table. 

In  this  paper  the  actual  procedures  for  matching 
segments  were  not  completely  addressed.  The  above 
examples  suffice  to  point  out  that  comparison  based  on 
curvature  of  segments,  or  other  space  invariant 
properties  of  the  segments  is  also  an  integral  part  of 
the  partial  shape  recognition  procedure. 

VII. b. ii  REASSIGNMENT  AND  RELOCATION  OF  CRITICAL  POINTS 

BASED  ON  SYMMETRY 


There  is  one  more  problem  associated  with  the  case  of 
the  circular  arc.  It  is  true  that  most  humans  assign  the 
highest  degree  of  membership  to  the  end  points  of  the 
arc  but  it  is  also  true  that  most  of  them  do  not  assign 
similar  degree  of  membership  to  the  remaining  points  on 
the  arc.  Specifically  it  is  usually  assumed  that  the 
critical  points  found  relative  to  the  axis  of  symmetry 
are  more  important  than  the  rest.  It  can  be  argued  that 
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the  determination  of  such  an  axis  of  symmetry  occurs  at 
some  post-cognitive  level.  In  other  words  it  occurs  only 
when  parts  of  the  shape  are  'compared'  with  itself.  The 
comparison  signified  by  the  word  'compared*  in  the 
previous  sentence  occurs  before  the  axis  of  symmetry  has 
been  located.  It  is  this  kind  of  comparison  which  has 
been  addressed  in  this  paper.  This  of course  excludes  the 
redefinition  of  critical  points  and  r assignment  of  their 
degrees  of  membership  based  on  the  feed-back  of  the 
results  of  comparison  of  the  shape  with  a  copy  of  it¬ 
self  (see  Fig.  1.),  and  falls  outside  the  dotted  line 
shown  in  Fig.  1. 

Vll.b.iii  HOW  TO  SET  THE  THRESHOLDS. 


An  important  question  associated  with  the  above 
method  is  that  of  the  thresholds  in  the  algorithm.  More 
features  can  be  determined  at  lower  thresholds,  while 
many  features  can  be  eliminated  at  higher  thresholds. 
Resampling  the  curve  at  uniform  arc  length  is  a 
necessary  first  step  in  the  determination  of  the 
thresholds;  even  though  we  did  not  resample  some  of  the 
shape  at  uniform  arc  length  to  get  a  feel  of  the 


robustness  of  the  algorithm.  This  is  the  reason  why 
some  of  the  critical  points  were  missed  near  the  leg  of 


the  pig  shape  in  Fig.  16.  In  general  the  exact 
threshold  will  depend  on 

1)  Resolution  of  the  Machine  (  a  machine  may  have 
variable  range  of  thresholds) 

2)  Knowledge  from  the  data  base  (see  Fig.  1.) 

3)  Feedback  from  the  output  table  (see  Fig.  1.) 

It  is  not  clear  how  the  thresholds  can  be  adjusted 
using  the  feedback,  however  there  appear  to  be  at  least 
two  factors  which  seem  to  influence  this  aspect,  namely, 

1)  Error  or  Noise 

2)  Number  of  'small  features' 


The  former  means  an  estimate  of  the  error  between  the 
best  fit  polynomial  and  the  linearly  interpolated  curve. 
The  latter  comes  from  a  general  observation  that  humans 
tend  to  ignore  (average  out)  features  if  they  are 
smaller  than  the  overall  global  size,  and  if  there  are  a 
large  number  of  them.  However  if  there  are  only  a 
couple  of  small  features  they  often  become  a  point  of 
'focus'.  Thus  small  features  cannot  be  generally 
ignored. 


VII. C  FUTURE  WORK 


It  should  be  very  clear  from  the  problem  statement  in 
the  introduction  that  only  the  problems  in  the  area 
enclosed  by  the  dotted  line  in  Fig.  1.  have  been 
addressed  in  this  dissertation.  A  lot  of  work  still 
needs  to  be  done  in  the  future  before  the  complete 
implementation  of  the  system  shown  in  Fig.  1.  becomes  a 
reality.  Some  of  the  important  topics  that  need  to 
worked  on  are  listed  as  follows: 

1)  Edges  of  objects  or  shapes ,  have  been  successfully 
detected  in  images.  For  a  brief  summary  of  the  methods 
used  to  detect  edges  the  reader  is  advised  to  refer  to 
[16].  The  problem  of  transforming  the  edges  into  a 
sequence  of  points  defining  a  shape  has  also  been 
attempted  [10] ,  [16] .  These  methods  work  reasonably 
well  on  ideal  data  specially  when  the  points  are  located 
at  nearly  uniform  intervals  (e.g.  eight  connected 
neighbourhood) .  However  these  methods  fail  to  work  on 
general  edge  data  (i.e.  data  defining  the  edge  of  the 
object) t  which  is  normally  obtained  using  edge  detectors 
[16] #  on  actual  'real  life'  images.  So  work  needs  to  be 
done  either  to  alter  these  methods  so  that  they  can 
adapt  to  the  edge  data,  or  the  methods  of  detecting 


edges  need  to  be  modified  so  as  to  produce  output  data 
which  is  close  to  ideal. 

2)  The  question  of  thresholds  which  has  already  been 
addressed  also  needs  to  be  worked  on.  Explict  standards 
defining  the  range  of  resolution  of  the  machine  need  to 
to  laid. 

3)  Host  of  the  state  of  the  art  work  [30]/  [54],  [68], 
[71],  [88],  pertaining  to  the  interpretation  process 
described  in  Chapter  I.  deals  with  the  extraction  of 
motion  information  from  changes  in  the  projection  of  the 
object  on  a  plane.  Practically  no  actual  work  has  been 
done  which  exploits  the  relationship  of  the  object  with 
respect  to  the  background  towards  this  end. 

4)  Though  some  guide  lines  can  be  drawn  from  the  field 
of  psychology  for  the  implementation  of  a  workable 
knowledge  data  base  shown  in  Fig.  1,  the  actual 
implementation  of  a  general  of  a  data  base  comparable  to 
the  human  vision  system  is  still  a  far  cry  form  reality. 
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APPENDIX  A 


DETAILS  OP  THE  PROGRAM  FOR  THE  ADAPTIVE  LINE  OP  SIGHT 

METHOD. 

Itefering  to  the  flow-chart  shown  In  PIG*  27.  side  A  rtfsrs 
to  ths  condition  whan  ths  computations  ars  performed  from  I 
to  J  modulo  ths  number  of  points  in  ths  shape,  while  side  B 
refers  to  the  condition  when  computations  are  performed  from 
J  to  to  I.  Forward  track  or  Ptraek  denotes  to  a  condition 
when  X  is  held  constant  while  J  is  incremented.  Backtrack 
or  Btrack  denotes  a  condition  when  J-l  is  held  constant 
while  X  is  decremented. 

Initially  X  and  J  are  always  Chosen  to  be  adjacent 
paints,  first  going  in  the  clockwise  direction  then  in  the 
anti-clockwise  direction  The  details  of  the  computational 
block,  the  detection  block/pcocess  and  some  other  blocks  are 
as  follows 
COMPUTATIONAL  BLOCX 

Find  the  equation  of  the  straight  line  L  joining  I  to  J 
Find  the  distance  DXSTAXJ  between  points  X  and  J. 

Find  the  equations  of  the  straight  lines  normal  to  L  and 
joining  every  point  P  in  between  X  and  J. 

Find  the  intersection  (XXNTSEXJ  ,  XXNTSEXJ)  of  each  of 
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the  normal  linos  found  in  the  above  step  with  the  straight 
line  L 

Find  the  distances  DISTXIXJ  and  SISTYXXJ  from  point 
X  to  the  intersection  and  from  point  J  to  the  intersection, 
respectively  for  every  P. 

Find  the  normal  distance  NDISTAU  from  every  9  to  the 
straight  line  L. 

rind  the  normal  vectors  from  the  straight  line  to  every 
point  9. 

DETECTION  PROCESS 

In  this  block  a  detect  switch  is  set  indicating  that  a  crit¬ 
ical  point  has  been  detected  at  X  and  J-l,  if  at  the  first 
instance,  a  point  is  found  which  is  not  on  the- same  side  of 
L  as  other  points,  or  a  point  cannot  be  mapped  injectively 
on  to  the  straight  line  L.  The  former  condition  is  checked 
by  comparing  the  magnitude  of  the  sum  of  every  two  adjacent 
normal  vectors  with  the  magnitude  of  those  forming  the  sum  + 
a  THRESHOLD!,  while  the  latter  condition  is  checked  by  com¬ 
paring  the  sum  of  the  distance  DXSTXXXJ+DXSTYXXJ  to  DISTAXJ 
+  THRESHOLD 2.  Where  thresholds  1  and  2  are  set  to  account 
for  round-off,  truncation,  quantization  and  other  errors. 
SWITCH  BOARD 

This  is  a  control  block  which  forces  the  computations  to  oc¬ 
cur  in  an  alternating  sequence  FORWARD-TRACE  -BACK-TRACE 


m  TBS  WO  RELATIONS? 

In  path  B  the  two  eolations  aeo  not  satisfied  if  ono  of 
tho  following  is  truo. 

1)  Tho  x  oe  y  eoapononts  of  tho  intersections,  i.e., 
XINTSEU  and  TINTS  CU  respectively,  are  not  no  no  tonic. 

2)  Tho  son  of  tho  distances  DISTXIU  +  DISTTIU  is  not 
equal  to  tho  distance  CISTAU  ±  THRES10LD2 


APPENDIX  B 


DETAILS  OF  THE  FEATURE  VECTOR  FORMATION  PROGRAM. 


To  find  whether  a  point  within  (n  +  3}  of  the  nth 
critical  point  belongs  to  the  nth  feature  it  is  not 

necessary  to  carry  out  the  computation  for  each  of  the 

points  within  (n  +  3) .  Specifically  the  following  two 
facts  have  been  exploited  to  reduce  computation: 

1)  Points  within  (n  +  1)  are  always  in  sight  of  Cn 

2)  The  Line  of  Sight  of  a  Point  is  a  symmetric  relation. 

Let  the  relation  (is  in  Line  of  Sight  of  the  Point) 

be  denoted  by  R;  then  if  C1  RCj,  implies  R  q  .  Once 

it  is  determined  that  R  Ch  then  it  is  not  necessary 
to  compute  if  R  C|,  In  context  withr  our  example, 

c  n  R  c  n+2  C„t2  R  Cn_2 

and 

C  n  R  C  n+3  implies  Cn+3  R  Cn_3 

Thus  it  is  necessary  to  check  only  if  the  critical 
point  Cn+2  and  Cn+3  are  in  Line  of  Sight  of  Cn. 

If,  the  point  of  intersection  of  the  line  joining  Cn 
to  Cn+k,  k«2,  3.,  with  a  line  between  Cj  and  C j+j 

j*l...N. r  is  denoted  by  INT,  and  the  coordinate  of  the 


following  points  are. 


INT  »  (XINT,  YINT) , 

Cn  -  (FEAT1X,  PEAT1Y) , 

Cn+k  -  (PEAT2X,  PEAT2Y) , 

Cj  -  (CDRVE1X,  CORVE1Y)  , 

Cj+1  “  (OTRVE2X'  CURVE2X)  , 

and  the  distances  between  the  following  points  are, 

INT  and  C  »  DIF1, 
n 

INT  and  Cn+k  -  DIF2, 

INT  and  Cj  ■  DIC1, 

INT  and  Cj+1  -  DIC2, 

CR  and  Cn+k  -  DISFEA, 

Cj  and  Cj+1  »  CURVED, 

then  the  relation  R  is  not  true  if 

DIP1  +  DIF2  -  DISFEA 
or 

DIC1  +  DIC2  «  CURVED. 

The  matrix  Stf SET  in  the  program  is  this  binary 

relation  table.  The  output  of  the  first  proc  matrix  has 

seven  columns.  The  first  column  corresponds  to  the 

sequence  number  of  the  critical  point  in  the  critical 

point  list.  The  next  two  columns  are  the  coordinates  of 

the  critical  points.  The  last  four  columns  contain  the 

boolean  relations  operator.  Specifically  they  hold  the 

following  the  binary  relations; 
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v.v.vncxv;-. 


APPENDIX  C 


DETAILS  OP  THE  PROGRAM  FOR  POST  PROCESSING. 


This  program  takes  as  input  the  critical  points  which 
are  represented  by  their  sequence  number  in  the  shape 
data  list,  and  process  them  sequentially  according  to 
the  following  steps: 

1)  Every  set  of  five  consecutive  points  is  replaced  by 
their  median  point. 

2)  Every  set  of  four  consecutive  points  is  replaced  by 
the  second  point  if  PRC  «  SET  otherwise  by  the  third 
point  if  FLC  -  SET.  Both  PRC  and  FLC  cannot  be  equal  to 
SET  at  the  same  time. 

3)  Every  set  of  three  consecutive  points  is  replaced  by 
their  median  point. 

4)  Every  set  of  two  consecutive  point  is  replaced  by  the 
first  point  if  TLC  ■  SET  otherwise  by  the  second  point 
if  TRC  -  SET.  Both  TLC  and  TRC  cannot  be  equal  to  SET 
at  the  same  time. 
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PROGRAM  LISTING  FOR  THE  ADAPTIVE  LINE  OF  SIGHT  METHOD. 

INPOT:  Three  column  data  set  named  SHAPEXY.  The  first 
column  represents  the  sequence  number  of  data  points  in 
the  shape  data  list.  The  2nd  and  the  3rd  columns  are  the 
x  and  y  coordinates  of  the  shape  data  points. 

OUTPUT:  1)  A  10  column  matrix  array  named  INF.  The 
columns  contain  the  following; 

1st  column  ■  The  number  of  critical  point  pair  (the  Ith 
and  the  Jth  point)  found  in  the  first  pass. 

2nd  column  «  sequence  number  of  the  Ith  critical  point 
in  the  shape  data  list. 

3rd  and  4th  columns  »  the  x  and  y  coordinates  of  the  Ith 
point  in  column  2. 

5th  column  «  sequence  number  of  the  Jth  critical  point 
in  the  shape  data  list. 

6th  and  7th  columns  ■  the  x  and  y  coordinates  of  the  Jth 
critical  point  in  column  5. 

8th  column  ■  the  location  (sequence  number  in  the  shape 
data  list)  of  a. maxima  or  a  minima  in-between  the  Ith 
and  the  Jth  point. 

9th  and  10th  columns  ■  the  x  and  y  coordinates  of  the 
point  whose  sequence  number  is  in  the  8th  column. 


OUTPUT  2)  The  input  data 


•••••••••START  SYNCHRONISE? 

PROC  MATRIX  ? 

FETCH  SHAPEXY; 

MH-MROW(SBAPEXY)  ? 

RESET-#; 

SET-1? 

FTRACX-SET ; 

BTRACX-RESET? 

DISCON- RES  BT ; 

CONNEC-SBT? 

PATHB— DISCOR ; 

SIDEB-SET; 

0PSXG8T— SET? 

**•*••*••••*•  STARTING  AND  ENDING  POINTS  FOR  THE  ALGORITHM  ARE  I  AND  J  ; 

I-ll 

J-2? 

ITEXP-I? 

JTEXP-J? 

NN— 1; 

••••••••••••«  NN  IS  A  INDEX  WHICH  POINTS  TO  THE  NEXT  CRITICAL  POINT; 

•••••••••••«•  NEED  A  MATRIX  TO  STORE  INFORMATION  ABOCT  A  CRITICAL  POINT 

■ay  tha  axpactad  matbac  of  critical  point*  i<  75  ao  diaanaion  INF; 
INF-J.  (73.11.0? 

•? 

•? 

•? 

•? 


•COMPOTAT  ION  BLOCX; 
START? 

DETECT1- RESET; 
DETECT2-RESET; 
DETECT3- RESET; 
DETECT-RESET? 
EXIT-RESET? 


AQ-SHAPEXY ( J.3)  -SHAPEXY  { I.  3)  ? 
BQ-SHAPEXY  ( J.2)  -SHAPEXY ( 1 . 2)  ? 
CQ-SHAPEXY ( J.2) ISBAPEXY (1.3) I 
DQ-SBAPEXY (1,2) PSBAPEXY (J.3) ? 
SLOPE  LJ-AQ#/BQ; 

INTER U-  (CQ-DQ)  f/BQ? 
DISTAL3-S0RT(  (BQtBQ)  +(AQ»AQ) )  ? 


MCE  2 


•MIIE  TO  KNOW  HOW  many  mints  ARE  IN  BETWEEN; 


•l 

if  (i>j) 

THEN 
DO  I 

PQ-NH-U+l)  > 
PT-J+PQ; 

END  | 

ELSE 

DO; 

rr-j-tifl) i 
END; 
ir  pt-b 
THEN 
CO  I 

GO  TO  SNBOAROi 
END; 


DIMENSION  THE  MATRICES  FOR  MINTS  IN  BETWEEN 


NINTERCT*J. (PT.l.B) 
XINTSELJ-J. (PT.l.B) 
YINTSBLJ-J. (PT.l.B) 
KDISTALW. (PT.l.B) 
CISTXIIJ-J. (PT.l.B) 
CISTYIIJ-J. (PT.l.B) 
VECTORXI-J. (PT.l.B) 
VECTORTI-J. (PT.l.B) 
SUNVECXI-J. (PT.l.B) 
SUMVECYJ-J. (PT.l.B) 
HAGVECXX-J. (PT.l.B) 


CERTAIN  OPERATIONS  HEED  TO  BE  PERFORMED  BETWEEN  MINTS  I  AND  J 


I 

CONSVEX-B ; 
CONSVEY-B  > 
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mt  3 


00  «•!  *0  ft  BY  1| 


XI-E+tfUl; 

ZZ«NW(ZZ«m)i 

it  (zz>«)  the*  do t  zz-ms  no* 

ASLOPE  ZJaMS  (SLOPE  IJ)  » 

ZP  (ASLOPEZJ  <  0  .0000001) 

nn 

VO, 

XZWSEZJ(X,)«SBAPBXYtZZ,2)  » 

YIMTSEU  ( K  r )  “SHAPEXY  ( I » 3 )  > 

EMD> 

ELSE 
00 1 

m— 10/slorzji 

lormcTd, )  ■hapezy(zz#3)  -(shapext (ii,2)  pm>  , 

XZMTSEZJ  (E, ) •(AIHTEECTIX, )  -ZNTXXZJ)  0/ (SLOPEZJ-IS)  > 
IZRUZJd, )  -(M0XZRSEZJ  (E, )  +NXXTEKY(E, ) )  » 
noi 
•i 
•» 

AQ«S«APEXY(II,2)-XZlTfEZJ(E,)  ; 

BQ-SlAnXY(ZZ,3)-YZirtSEU(E>)  > 

CQaSHAPEXY ( I , 2 ) -XZMTSEZJ (E, ) I 
DQ-S1APEXY t I . 3 ) -YZ1TBEXJ ( E , ) , 

PQ-nAPSXY(J,2) -XZHTBEZJ  (E, )  » 

G0-E1APCXY  ( J,  3 )  -YZHTSEZJ  ( E, )  , 

*» 

*> 

DZETXIZJ ( X. ) -1Q1T ( { CQ 0CQ) ♦ ( OQ 0DQ) ) , 

DXSTYZZJ  (E, )  »SQRT(  (FQ0FQ)  +  (GO»GQ) )  » 
■DXSTAXJ(E,)«SORT(AQ0AQ*BQ0BQ)  , 

VXCTORXZ (X, ) »AC; 

VECIOEYZ ( E , ) »BQ j 

SUMVECXKX,)  -COHSVEXt-VECTOAXZ  (E, }  ; 

SOWVECTJ  (E, )  -COKSVEY*VECTO*n  (E, )  ; 

AAGVSCXYtE, )  •SO«T(  (EOIWECXI  (X, )  0*2)  ♦(SOHVECXI  (E. )  *02) )  » 
COHEVEX-VBCTORXKE,)  * 

OONSVEY-VECTORYZ (E, ) » 

•l 

*, 

*1 

•» 

*1 

•» 

EXOl 

•l 

*» 

•> 
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it  (MnsaoucoN)  nn  go  to  sxoea; 
it  (sxon-sn)  nn 

DO) 

lire  axartr 

xr  (omxgbt-reset)  nn  go  to  quit; 
xr  <onxGit-«R)  nn 

DO; 

Un  LOCATES; 

exit- jet > 

GO  to  QOXT; 

EHOt 

no> 

3  IDE  At 


orra  detect  processes 


•••••••••••••OIVBKS  tO  DETECT  PROCESSES  1  AMD  2  •**•••*•••••••*•*•• 

•l 

*» 

•  PROCESS  DKTECT1 > 

DCnCTX  ■RESET; 

DO  14  10  R  IT  1; 

LL-DISTAXJ-».fJ; 

CL-DISTALJ  ♦».»»; 

CL-DISniXJ  (I, )  ♦DISTIIU  (E, ) ; 

IP  ( (CL  <  LL>  OR  (CL  >  OL) ) 

nn 

DO; 

DETECT1-KT; 

GO  TO  CONVERGE ; 

CRD; 

•; 

•; 

•; 

END; 


•; 

*; 

•; 

•  DETECT  PROCESS  2; 

*; 

COMMORXT-I.S; 

•; 

•; 

00  *»1  TO  Pt  ay  i; 
AQ-MDXSTAXJ(E.)-2>; 
BQ*COItMORXT  *29; 
COIWORET-KDXSTALKK. ) ; 


mi  s 


IF  ( (NKVtCXT(K, )  <  AO)  OK  (HASVECXTtK. )  <  KQ) ) 

Tan 

OOf 

DETBCT2-SET) 

60  TO  COHVIKll 

no> 

•» 

•» 

•> 

UDi 


•••» 

•» 

*> 

•j 

•i 

cowmeit 

*» 

•» 

DETECT- BESET; 

ir  {(omen  eq  art)  ot  tPKTtera  eq  set) )  Tan  dstect-set; 
IF  (DETECT  KE  SET) 

Tan 

DO) 

60  TO  SmOMtDt 
BED) 

rr«» 

OOf 

IF  (FTKACK-SET) 

Tan 

OOf 

oujrr-FT; 

AX-I» 

■I-J-1HW) 

BZ-NOO(BX,HM)  ; 

IF  (BZ-S)  Tan  BZ-NN) 

OLDI-AZf 
LINK  LOCATE I ) 

UD  LOCATE!) 

LHU  STOKE  I ; 

•Lin  cases  i 

IP  (EXIT-SET)  Tan  60  TO  QOIT; 

60  TO  SMSOMO) 

IF  (BTKACX-SET) 

Tan 

DO) 

antr-FT) 

M«NOD(Al,m)  ) 

IF  (ASH)  Tan  AZ-IM) 

anx-Mi 

BJ-J; 

IF  (OLOI  EQ  NEW I)  THEN  60  TO  SWBOAKO) 


LOX  LOCATZI; 

Lon  LOCATXX) 

LUK  JTORII; 

ir  (ixiT»arr)  nn  so  to  oorti 
un  DCT3  j 
IT  WttCT3»«*T 

n 

00 1 

im  c i*a. 1 1 

SO  TO  (HBOABOr 

no  I 

ir  ornen-mn 

nn 

DO) 

UB  ALTAI ; 

00  TO  AIOAIO; 

no  I 

OB) 


VKTOOXX 


nratu 

■unu 

oiraiu 

OUTTZU 


> 

•  ntm  boausi 

•••••••••«••••••••••■•••••••••••••••••••••••••••••••••••••••••••••••• 

SNOOMOi 

•» 

*» 

•l 

ir  t  (Drier  to  user)  *  irmac  ao  ird  - 
nn 
00 1 

J-MOO(J,MR)> 

if  <jm)  nn  j-m> 


VAC  I  ? 


xr  ( (DCnCT-KBID  t  (STHACS-SIT)  > 

TIM 
DO  i 

I- 1-14401) 

X-H0D(X,IUI)  ; 

xr  (X-S)  BO  x-m> 

CO  to  STMT) 

no) 

•t 

*; 

•) 

xr  ( tPBTBCT  BQ  SR)  A  (RXACX  V)  SR)) 

tin 

00) 

rBTMCS>SR) 

mAOMtSSR) 

BTIACS-SR) 

J-J-l'HOI) 

J-MOO(J,BM)  ) 
xr  (J-i)  Tin  J-M) 

X«fl<MM) 

X-MOD(X.M)) 

xr  (i-s)  tin  x-ni 

00  TO  STMT) 

no> 

xr  ((OSTSCT-SR)  t  (BTMCS-SR)) 

tin 

DO) 

SmtACK-SR) 

riacs-sR) 

STBACX-KSSR) 

I-J+W) 

x-Bood.m  i 

xr  u-i)  tan  x-m; 

j-14-i+nt) 

J-NOO(J.HM)  ) 
xr  (j-«)  Tin  j-wi) 

xr  (pathb-oiscoi)  tan  co  to  stdaxcr; 
xr  (sxon-usR)  tin 

DO; 

sxon-BR) 

XTECP-X) 

JTSXf-J) 

I'JTDi 

J'XTU) 

BID) 

00  TO  STMT; 

QSTOOTi 

SXDn-MSR; 

i-rm»; 

JaJTSXP; 


that  • 

sntAiesTi 

60  tO  (TAM; 

CRD» 

I 

•  un  rootxmcc 


LOCATtli 

DO) 

A-MDISTAIJ  (1, )  ) 

LOC-1) 

MQ-PT-1) 

DO  III-2  TO  RQ  >1  1) 
zr  (A  <  RDISTAU  ( III, )  ) 

na 

DO) 

A-MOIXTAU  (III, )  ; 

LOC-IIIl 

CRD) 

CRD) 

RMORM; 

•) 

•) 

•» 

LOCATE! i 

ir  ((M<4)  OR  (RR>7»))  THEM  60  TO  JOHTS) 
•AX-I+l+l M) 

•AX-HOD (AX, HN) ; 

•IP  (AZ-C)  RO  AX-NHi 

•BZ-J-l-Hm)*BZ-flOO(CX,mi)  f»IP  (•!«•)  RB  XZ-NH; 
MQ-PT-1  lHQO-PT-2; 

MAM-RDISTAXJ(liMQ,)  ; 

MRN-RDXSTAIJ  ( 2 1 PT, )  ; 

DAAD-RAM-iaH; 

DRAD-0ABD(2tNQ,) » 

DRAO-OBADPDAAD  ( 1  <  HQO . )  ) 

DAASD- A1S ( DEAD ) ) 

COOMTNX-4; 

DO  IXMI-1  10  MOO) 

AJOU-ORADtlXXXI.)  | 

CDIE-OAASO(XXXXI,)  ; 

IP  ((AHU<«)  AMD  (SKMS>4))  TMEH 
DO) 

IMP (MM, 1) •MR; 

DIP  (HR, 2)  -AX; 

IMP  (MM,  3  i  4)  -SHAPEXT  (AZ  ,  2  1 3)  I 
IMP  (HR,  S)  MX) 

IMP  (MM,  4  1 7)  -SHAPEXT  ( BX ,  2 1 3)  ; 

COOMTNI-COOMTNE+1 ; 

HtLOC- KXXI+NM*  AX  ^l) 

KTLOC-MOD (HtLOC, MR) i 
IP  (HTLOCM)  THEM  MTLOC-MM) 

•IP  ((AIM  <  I)  AND  (SMI  >  .27))  TIXM  DO) 

IMP (HR, •) -HTLOC; 


MSB  » 


ar(m,»u*)>aumx(imoc,2si)  > 

i r  (M»7f ) 

nn 

00 1 

BUT*  SRl 

60  TO  JOIM* 

DO) 

noi 

■mi 

JOHN) 

tr  <m>7«)  nn  «xit-sbt» 

BROW) 

*) 

•) 

*) 

•) 

STOBBIi 

mr(H,i)-n> 

xar(w,2i-U) 

W(W,3i«)-SBAr«T(A2,2t3)  > 

nr(w.S)-U) 

XBP(W.«i7)  *CBABSXY(BZ>2i3)  ) 
LOOLOOBH+AI  l 
LOO-MOD(bOC.M)  ) 

ir  (loc  •«)  nn  loc-w> 
nrtm.D-LOC) 

IBP(W.»«X*)-«BB»*n(tOC,2<3)  ; 
ir  (W  >  7«  ) 

nn 

00) 

niMiii 

no) 

BB-BB+1) 

BROW) 

•) 

•) 

•) 

CIRCLE t 
BROW) 

•) 

•) 

•) 

CBECBj 

BT-BB-4) 

00  BX*  1  TO  W) 

AJt-nr(Kz,S) ) 

BBBZ-BX+2) 

BBBY-BX-2 ) 

rr  (BBBT  <  AX  <  BWZ) 

nn 

DO) 

CUT*  BUR) 


PASS  1 • 


BMD) 

inwi 

•» 

•i 

•> 

•XUTl 

onxsn-UTi 

•If  you  mum  s*ro  ehraahold  chan  uaa  tha  go  balov 
othar  viaa  daiata  tha  go  to  itAtaoant. ; 

eo  to  xas> 
orsxen-usttf 

ip  (  (pt>«)  oi  (!■>?■) )  tan  oo  to  jokpm » 

nm  DIITAXJ  l 
DO  XB-1  TO  FT  IT  1 1 

u-oianoj-.xi 

ouDxsmxj+.xi 

CL-Dirrxilj  («, )  ♦  0ZCT7ZIJ  (TI, )  I 

non  oi  u  a.i 
ir  aa<ui  ob  (a.>oi) ) 

tan 

00 1 

■otaro-nt; 

caoi 

KLSBi 
DO  I 

0PaiQBT-8BTl 

00  TO  JOMPH; 

no  i 

aao  i 

Hfi 

MQ-FT-l>MQQ-PT-2» 

TXAT-xirraiu (x>mq,)  i  t2ax»xikt»iu<2xPT,) » 

TAX-TIAT-T2AX ; 

TAlX-XAXUiKQQ,) IIA2I-XAI ( 2  tKQ, ) t 

zaax-xaxx*xax x» 

xxax>xxataaiJ(XiaQ,) »i2bi«iibtbxij(2ipt,) » 
xax-xxax-xaaxi 

xaxx-xaxdiaoo,)  ixa2x-nx(3iNQ,)  r 
tbit-tbit  mix  t 
nan  xaax  taax  » 

DO  XHB-1  to  MQQ; 

xtxoo-xaax (die,  )  iXLXoaaxiBTinu, )  i 

ip  (  (xkxoo  <  .ttgggg)  oi  (xixoo  <  .•••■••>)  tan 

DO  I 

optiorr-Aurri 
no  i 
pro  i 
JOKPM t 

anon  * 


ALTtJUs 

m-HH-ll 

M3-m*2> 

*-  nr  («i,  5)  * 
»-ur(«2,S) , 
Tnin-mr  («i.) ; 
IP  (A-tl 
Tin 
DO  I 

nroa.i-nwii 

!W-SX ; 

am* 

unni 

•t 

•; 

*; 

DCT3 1 

izzi»up(w,2)  i 
izxj-uF(mi,3)  > 
vnm-m-ii 

ir  azzz  <  tzzj) 


PAG  i  u 


oo  us>i  10  aomi  sx  i» 

IF  ((IZZX  <UF(IZZ>2)  <XZZJ)  01  (XZZJ  <IKP(IZZ,S)  <IZZJ) ) 

no 

ro» 

ocncz3-m» 

own-xzzi 

Ft 1ST  DKTCCT3  OHOTt; 

not 

•» 

ns) 

•l 

IP  (XZZX  >IZZJ) 


00  IZZ-1  TO  MUM  IT  1; 

UDCXX*  UP  ( IZZ  >  2!  ; 

X*D«XJ-UP(IXX,J>  I 

IP  (((XZZX  <  XHDtXI)  Ot  (XMDCXX  <XZZJ) )  Ot 

((xzzi  <  itozxj)  ot  (xzzx  <  iHoexjm 

nn 

DO; 

0ZTtCT3«»fT> 

0«0TX-IZ1) 

FXXRT  0CTCCT3  IZZ; 

C8D; 

(VO; 

U101UI; 

QOITi 

OOTFOT  8SAFCX7  OUT-SHAPCXT  (*niArt*-  (C0L1*J  C0L2-X  COU-X) ) ; 
OOTFCT  IMP  OCT» UP ( tCKAXta ( C0X>1  •  J  C0L2-II  C0L3-XII  C0L4-XII 
C0L5-JJ  C0L4-XJJ  C0L7-WJ  C0L8-L0C  COL 9- XL OC  COUI-TLOC) )  ; 


FEATURE  VECTOR  FORMATION  PROGRAM 


INPUT:  The  input  to  the  program  is  data  set  called 
SHAPEXZ .  SHAPEXZ  contains  the  X  and  Y  coordinates  of 
the  critical  pc  nts. 

OUTPUT:  The  output  of  the  the  program  is  a  matrix  called 

DICTION.  DICTION  has  seven  columns.  The  1st  column  is 

the  sequence  number  F  to  which  the  critical  point 

n 

belongs.  The  2nd  column  is  the  sequence  number  of  the 

critical  point  in  the  list  of  critical  points.  The  3rd 

column  represent  the  binary  relationship  C  R  C  ,  k»l, 

n  n+k 

2,  3,.  The  4th  and  the  5th  represent  the  x  and  y 

coordinates  of  the  critical  point  if  R  «  1;  The  6th  and 
the  7th  column  contain  the  feature  vector  components 
and  a  nk  respectively. 


PHOC  NAT* IX  ; PITCH  SHAFIXZ  ;H*-N*ON( SHAPtXZ)  ; 

»■«) 

COH-J.  (HN.7,.) ; 
mu-nt-ii 
nbpi-nh-i  i 
nwu-nn-1; 

COH  <  1  ill,  1 1 2)  •SHAPtXZ ; 

C0N(BHP1;NN,1;2) -SHAPtXZ (1 i «,) ; 
COB(liNN>ll,3i«)-CO*(2;NN,lt2> « 

COH  (Ml,  3 1 4)  -COM  (7,1 1 2)  ; 

C0N(1 3HH,  5  i 6)  — C3NI1  ;NH,  X  x 2 )  -CON ( 1 :NN, 3 :4)  ; 

G0N(liNN,7)  -CON ( 1 iNN, S )  #CON( liNN, 5)  ♦CON  ( 1  :NN,C)  ICOHd  >NH,() 
ABIU-CON(  1  zHM,  7)  > 

COH(llNN,  7)  -SOKT(ABNE)  ; 

•» 

•» 

*1 

DO  HA  -1  TO  KM  i 
AB CD-CON (HA, 3 ) , 

ABCD— ABS(ABCD)  ; 

IP  (ABCD  <  .IIHMl)  thin 
DO  I 

COR(HArS) -9999999; 

COH (HA, 6) -SSISSS; 

nti 

IP  (ABCD  >  .■•■••11)  THEN 
DO  ; 

CON(HA,3)-CON(HA.<)*/CON(BA,3> ; 

COH (  HA, 6) -CONI  HA,2)-COH(  HA, 5) ICON (  HA,1); 

END) 

END) 


NNM3-NN-3 ; 

NNH2-NN-2 ; 

mms-m-Si 

SHSCT-J.  (NN,  4  , . )  ; 

DO  1-1  TO  NN; 

IP1-I-1; 

IP2-X-2; 

IP3-I+3; 

•; 

•l 

COO  NT-1; 

RENTER i 

IP  (COO NT-1)  THEN  IP-IP2; 

IP  (COO NT-2)  THIN  IP-IP3; 

•> 

IP  (  IP  >  NN)  THEN  SO  TO  EXIT; 
PEAT1X-C0N (1,1) ;PEAT1Y-C0N (1,2) ; 


PACE  2 


PEAT2X-CQH(IP,1)  ;P*AT2Y-CON(IP,2)  ; 

•t 

*1 

P-CONU,)  ;G-COE(XP, )  I 
•PEUIT  P  G; 

•SQOATIOH t 

0ELPEX-PEAT1X-PEAT2X | DELPEY-PEAT1 Y-PEAT2Y ; 
DXSPEA-DELPEX 90ELPEX  ♦  DEL PE Y 9DELPE Y  j 
DISPXA-SQHT(DISPXA) ; 

DTEST-AES(DBLPEX)  ; 

IP  (WEST  •<  .MMM1)  THEN 
00; 

PH-99999999; 

rc-«; 

EES; 

ir  (OTEST  >  .•••(Ml)  TEES 
00; 

IW*0ELPEY9/0ELPEX ; 

PC-PEAT1Y-PH9PEAT1X ; 

EEC; 

•; 

•; 

•; 

•; 

SICHT-i! 

•smart 

00  BI-  1  W  EE; 

COEVEN-COB (BI. 3) 1 
COEVEC— C0E(EX,6)  ; 

COEVED— COE(EX,  7) ; 

Dtrsto-(PN-COEVEN) ; 

DIPSLOP— ABS  (DXPSLO)  ( 

IP  (DIPSLOP  >  .IIIIH1)  THEE 
DO; 

IP  (COEVEN  >-  9999999)  THEN 
DO; 

XIET-CON (NIil) ; 

YXNT-PH9XIHT*PCj 
GO  TO  STEPX; 

EMC; 

*; 

IP  (PH  >•  9999999)  THEE 
00; 

XIET-PEAT1X; 

YIET-COHVEN»XIHT*aiKVEC ; 

GO  TO  STEPX; 

EES; 

XIHT-(COEVEC-PC) 9/OIPSLOf 

YINT-PH9XXHT-PCf 

STEPX; 

CDHVE1X-C0N ( HI , 1 ) ;C0EVE1Y-C0N ( HI, 2) ; 
C0HVE2X-C0N(NX,3)  ;CUKVE2Y-C0N  ( NX, «)  ; 
PSTXD1-PEAT1X-CURVE1X; 


iii] 


FACE  3 


I 

FSTYD1-PEAT1 Y-CDRVE1 Y ; 

FSTXD2-PEAT2X-CORVE1X; 

FSTYD2*FEAT2Y-CURVE1Y  ; 

SECXD1-FEAT1X-CURVE2X ; 

SECYD1 »PEAT1 Y-CDRVE2 Y ; 

S8CXD2-FEAT2X-CORVE2X ; 

SECYD2-FEAT2Y-CURVE2Y ; 

*; 

*; 

•> 

I F  (((FSTXD1-0)  A  (FSTYD1-0)  )OR(  (FSTXD2-0)  MFSTYD2-0)  )OR 
( (SECXD1-0) i (SECYD1-0) )  OR  ( (SECXD2-0 ) A (SECYD2-0) ) )  THEN  GO  TO  OUT9 ; 

DIFXX-FEATIX-XINT; 

DIF1X-FEAT1 Y-YINT ; 

Diri-DIFIX  »0IF1X  +DIF1 Y »DI F1Y ; 

Diri-SQRT(DIFX) I 
OIF2X-FEAT2X-XINT; 

0IF2Y-FEAT2 Y-YINT ; 

DIF2"DIF2XtDIF2X+DIF2Y»DIF2Y; 

OIF2-SQRT<DIF2) ; 

DIC1X-CDRVE1X-XINT; 

DIC1Y-CORVE1Y-YINT; 

DIC1-DIC1X9DIC1X+DIC1YIDIC1Y; 

DICl-SQRT(DICl) ; 

0IC2X«CBKVE2X-XINTf 
DIC2Y-C0SWE2Y-YINT; 

DIC2-DIC2X#DIC2XfDIC2Y«DIC2Y; 

DIC2-SQRT(DIC2) ; 

OTHDIODIC1*DIC2;SOKOIF»DIFX*DIF2; 

OIFFl-aiRVED-SOfU3ICjOrFFX«A8S{OXFFX)  ; 

DIFF2-DISFEA-SUHDIF;  DIFF2-ABS (DIFF2) ; 

*» 

O-CONtNI.) ; 

•PRINT  0  XXNT  YXNT  DXFX  OIF2  DIC1  DIC2  DXSFEA; 

•  • 


•PRINT  0; 
•: 


•> 

IP  ( (DIFF1  <  .090911)  t  (DIFF2  <  .099001))  THEN 
00; 

SIGHT-9; 

APTEST-CONd.l  :2)  ; 

CPTEST-C0N(IP,1:2)  ; 

BFTEST-CONfNI,) ; 

•PRINT  APTEST  BPTEST  CPTEST;  ' 

GO  TO  0CT1; 

END; 

END; 

OUT! ; 

END; 

OUT! : 


A 
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•IP  ( ( (I-l)OR(I-3)OR(I-4)OR(I-31)OR(I-32)OR(I-33)OR(I-34)OR(I-3f) ) ) 


•00 1 

•PRINT  APT  1ST  BITEST  CPTEST  0  XIHT  TINT  DIP1  Dir 2  Did  DIC2  DISFEA; 
•END; 

•PRINT  SIGHT; 


IF  tCOCNT-1)  THEN 
DO; 

SNSET(I,3)-3IGHT; 

COUNT-COONT+1 ; 

IF  (I  •<  MJIM3)  THEN  SWSET( IF2  >2)  "SIGHT; 
GO  TO  RENTER; 

END; 


IF  (COUNT-2)  THEN 
DO; 

3WSET(I,4) -SIGHT; 

IF  II  -<  mm3)  THEN  SNSET(IF3,1> -SIGHT; 
60  TO  LOO  POUT ; 

END; 


LOO  POUT: 

END; 

•> 

EXIT: 

•PRINT  SWSET; 


CON{4iNNM3,3iS)-SHSET(4iNNN3,) ; 
COM( It  3  >3 1<) -SWSET (NNH3 iNNN3 , ) ; 
CON(NNM2;RN,3tC) -SWSET(«:4,) ; 
FREE  SWSET; 

•PRINT  CON; 


OUTPUT  CON  OOT-CON (RENAME- (COU-II  COL2-I2  COL3-I3  COL*- 14  COLS -15 
COM-14  COL7-I7) )  ; 

PROC  PRINT  DATA-CON; 

DATA  CONI ; SET  CON DROP  17; 

PROC  DELETE  OATA-CON; 

PROC  PRINT  DATA- CO MI ; 

*; 

*; 


PROC  MATRIX;  FETCH  CONI;  NN-NROW (CONI ) ; 
NN7-74NN; 

NNNS-NN-S ; 


MSI 

coil  (ms  lira,  7)  -1/2/3/ 4/S/4; 

diction-j.  (mm,  7,.) » 
nou-m-i; 
do  at-  i  to  mil 
DTPl-XM+ls 

arn-ar-2; 
wrs-iwo  > 

DfM-W*4, 

DfPS-W+Ss 

nrx-i n*«; 

DfP7-IW*7 ; 

W7-W#7j 

nr7N4-nr7-4s 

W7IO-W7-3» 

Df7N4»W7-4s 

D»7K3-Df7-3; 

ar7«-nr7-2; 

nr7mi-nr7-i; 

xr  ((xwtc  >  mi)  ok  (inm  >  nn)  )  Tin  so  to  racks 
AWC-CON1 ( IWP3 ,3(4) ;  AW  r-COMl ( XWP3 , 5 1  < ) t 
DICTION ( IM7N4 i IM7 ,4 iS)  -C0R1  ( XW i  IWPt  ,1  ■  2)  S 
DXCTIOK ( IW7M I XW7,2) -CONI (XHi XWPC,7)  » 
DXCTXOM(IW7MiXN7,l)-XWP3«J.  (7,1,1)  s 
OICTXOV ( IW7K4  < IW7K2 , 3 ) — J . (3,1,1); 

DICTION  ( IW7M  tXW7KS,3)  -AMI'  s 

DicTioK(nnnitxin,3)-«wr'  > 

an; 

RACKS 

FRIZ  CO  Mis 

•s 

*» 

•s 

DO  Ul-1  TO  HMHs 
AFDT-DICTIOM (111,3) s 

xr  (AFDT-I)  TRIM  DICTIONt 111,4 s5> -J. (1,2,1) s 
ENDS 

PRINT  DICTXOHs 

•s 

•s 


•s 

•  SHAPE  VICTOR  CALCULATIONS  s 

•s 

*ss 

•s 

•s 

DO  1 33 1 -TO  HNs 
13317*133 

*S 
*S 
•( 

•; 


MCI  ( 


BO  X33X-1  SO  ■> 

X13X7-74X33X  »  X33XM-I31X7-4 1 

D-J.  (7,7, . )  ;XI-J.  (7,1, . )  ;TI«J.  (7,1, . )  ,XD-J. (7,1,.) ,DI«J. (7,1,.) ,PIT-J 
#1#*) > 

D-OXCTXON(X33XMi  133X7,)  ;XI-D(,4)  ,YI-D(,S)  ,M0-0(,3)  I 
XWNAX-XK+,)  »iXOIttY-YI(*.)  ,N-N0(+,>  » 

cent*oix«eionax4/n, 

CSMTBOX Y-SXSMAY  t/B  ; 

DX-XXtXX+YXIYX, 

OXCrXOB( X33XBB 1X33X7 ,f ) -OX i 

CAMC  LX- CENTRO  IXI/CXNTROIY; 

xr  (CXNTXOIX  >  •)  1BKB  RXr-AXCTAN(CANGLE)  , 

ir  tcxBWoix  <  •)  nxa  rxx-3.i«+a*ctaii(camgle)  » 

xr  ((  CEXTKOXX-X)  ABO  (CENTtOIY  >  «) )  THEM  KKT-1.57} 

xr  ( (CXNTXOIX-t)  ABO  (CEMTXOXY  <  •)  )  THEN  RET -4.71; 

*» 

•» 

•» 

*1 

00  IB  0-1  TO  7  » 

XB-X33XM+XB0-1, 
xr  (oxcsxaa(XB,3)-B)  nn 

DOi 

0ICSX0B(U,«)-#» 

XX(XBO.)-B| 

ITdBO.)-*; 

DICTION  ( IX,  7) > 

BHDs 

xr  (DICTION (IX, 3)  BE  •)  TIEN 
DO, 

fBT ( XBO , ) -XT ( XBO , ) #/IX ( IBO , ) , 

A- MY (1X0,)  , 

•t 

•» 

IT  (IX (IBO, )  >  •)  TIEN  DICTION (IX, 7) -AXCTAH (A) ; 
ir  (XX(IX0,)  <  I)  THEN  OXCTIONdH,  7) -3 .14*ARCTAN(A)  ; 

IT  ( (XX (IBO, ) — •)  AMO  ( IY ( XBO , )  >  ■)  )  THEN  DICTION (IH, 7) -1.57; 

IT  ((XXdBO,)-#)  AND  (lYdBO, )  <  *))  THEN  DICTIONdH,  7) -4. 71, 

DICTION (IH, 7) -REF-DICTION (IH, 7) , 

ACOXr-OXCTION (XX , 7) , 

Acoxr-cos(ACDir) » 

ASDXr-AXCOB(ACOir)  , 

DICTION (IX, 7 )-ASDir, 

END, 

END; 

non  o, 


PROGRAM  LISTING  FOR  THE  POST  PROCESSING  STEP 


INPOT:  A  data  array  named  AINF  which  consists  of  one 
column  of  numbers  representing  the  sequence  number  of 
critical  points  in  the  SHAPE  DATA  LIST. 

OUTPUT:  The  output  is  a  post  processed  data  array  AINF. 


For  more  details  see  Appendix  C 


HOC  NATH  IX)  FETCH  AIHF  iNAtXF'DkOW  ( AXHF)  ; 

HSSET-«)SET-1> 

FtC-SETlPLC-tESET) 

TlC-«tSST)TlC-SET> 

WTCOOWEl 

naimfs-maimp-i ;naixp2-naihp-2  )Naxmpi-mazmf-i » 
unr«*Mnr-«i 
•» 

so  mtn-i  TO  NAXMF4) 

HVW-MXV+4) 

ASNSSTOS-AXMPtMtWlHVV,) > 

AY001-AEME8T05(lr)  I AT002-AMISTO5  ( 2 , )  )AY003-ASXESTOS  (3, ) » 
AT004-A*KEST03<4,) >AT0CJ-A*MEST0S(5,) ) 
DTUl-ATOOl-ATOUaiOTUl-AM(DTDl)  , 

DTU2-AX002~AX003 | DTU2-AMS (DT02)  )DTO3-AI003-ATO04 ;DYU3-A1S(DTO3)  I 
DTO4-AI0O4-Ajr0OS)DT04-AMS(DTD4)  > 

ir  ( (0T01-1)  AMO  (DT02-1)  AMD  (DW3-1)  AMO  (DTD 4-1))  THEM 
00 1 

IHNP9-N1N+3 » 

A  IMF  (MW, )  -AT0C3 1 

AZMf  (HWU IMAXMF4, )  ■UEF (MWPS  iMAXMP, )  > 
AiMr(MAnniMAnr>)-j.(4.i,.) > 

(MO; 

KMDi 

•» 

00  M1M-1  TO  MAXMP3) 

WW-KTW-O) 

AEHEST04-AIHF(MWiKV¥,]  ; 

AXXX1-AKMEST04  ( it )  I 
AXXX2— ASHEST04  { 2, )  ) 

AXXX3-AHIUT04  ( 3  # )  I 
AXXX4— AEMEST04  ( 4 , )  I 

ADXXX1-AXXX1-AXXX2 ; ADXXX2-AXXX2-AXXX3 ; AOXXX3-AXXX3- AXXX4 ; 

A0XXX1- AMS  ( A030CX1 )  I A0XXX2-ABS  ( A0XXX2)  )AOXXX3-ABS(AOXXX3)  ; 
xr  ( (AD XXXI  EQ  1)  AMD  (A0XXX3  EQ  1)  AMO  (A0XXX2  BQ  1))  -THEM 
DO) 

MTVf  l-MW+l  )HWP2-MTVf2  ;  »YWP3-MW-t-3  )  MWP4-MW+4  ; 

XP  (FLC-SET)  THEM 
DO) 

AOIf  (MTM  * )  -AXXX2 ) 

EMC) 

IP  (PMC-SET)  THEM 
00) 

AXMP(NTN,)>ASXX3| 

EMD; 

AIMF  (HXVF1 :  NAIMP3 , )  -AXNP  (NTMF4  iMAIMF, )  | 
axmknaimf2«maimp,>-j.  <3,i,.)  i 


00  90  NAIMF2; 

HVY-HYMfl  | 

AXHEST02*AINP  (MW  iMVY, )  1 
ADDY-AXH8ST02  (1 , )  -AKH8S902  ( 2 , )  t 
AODY-ABS(AODY) , 

I F  (A00Y  BQ  1)  98  B) 

DO  I 

IF  (980889)  THE* 

DO  | 

AIMF(NYM,)-AHt88902(2,)  t 
ant 

HYNF1-MYW4-1  »MYWF2»*YW+2  > 

AZMPfHYWFl  tNAXMFl, )  »AISF(BY¥P2  iNAIMF, )  ; 
AXXFtRAXMF, )  aj.  (1,1,.)  I 

EMDj 

on  l 


009*09  A  IMF  009-AX*F(RDIAIU-(C0U«JJ)>  » 


PROGRAM  LISTING  OF  THE  REVISED  FEATURE  VECTOR. 


INPUT:  Data  set  named  AAUU  with  three  columns.  The  1st 
column  represents  the  sequence  number  of  the  critical 
points  in  the  critical  point  list,  while  columns  2  and 
3  represent  the  x  and  y  coordinates  of  the  critical 
points. 


OUTPUT:  The  program  has  two  output  matrices,  1)  STORES, 
2)  SKNO 

STOREE  has  six  columns  which  contain  the  following: 

1st  *  feature  number  F  to  which  the  critical  point  in 

n 

the  second  column  belongs. 

2nd  *  sequence  number  of  the  critical  point  in  the 
critical  point  list. 

3rd  and  4th  *  the  x  and  y  coordiantes  of  the  critical 
point. 

5th  and  the  6th  ■  the  distance  and  the  angle  component 

of  the  feature  vector,  namely  s  and  a  , . 

nk  nk 

SKNO  has  six  columns  which  contain  the  following: 

1st  ■  x  coordinate  of  the  mean  of  the  feature. 

2nd  ■  y  coordinate  of  the  mean  of  the  feature. 

3rd  *  the  distance  of  the  centroid  of  the  feature  from  a 


moc  MATS IX ; FETCH  AACO  ;*X-ltXO*  ( AACtJ)  JMH3-34NX-*  jSTORXE-J.  (NH3,S , . ) 
AAOO-5MAAD0; 

MO-mt-2; 

*» 

«> 

DO  1-1  TO  MOt 

in-i*iiir  (if i  >  ho  )  nn  im-im-mo; 

XH3-34Xt 
XH3N2-XN3-2; 
in*  1*2 1 

«TOM*(IM3M2tXM,l>-I»l»J.  (1,1.1)  I 
STOMB(XN3M2 i XH3.2 1 4) -AAOO ( I |XF2 ,1 13)  t 
XX-BTOXSB(XH3M2tXN3,l)  i 
TI-BIOKBB ( XM3M2 i XK3 , 4 )  1 
iam>a(M » 

IOIITI-YX(*,)  « 

HtMn-soim*/3> 

HtMnr-«nm«/i> 

CXMT-MAMMXAM  ♦  MAMYtNBMXt 
CBMTtO-HQMT(CBHT)  I 
XX-XK-MXANXM.  (3,1,1)  l 
TK-n-MAHTM.  (3,1,1)  I 
OX -XX*  XX  ♦  YItTXl 
DC-MQBT(DB) I 
*IXt-OB(f,>  t 

•TOftU  ( XH3M2 1 IH3 ,  S)  -DB#/SIXB  t 

•t 

•t 

It  (  (RXMX  -  •)  MO  (MAMY  >  •})  THEM  DO* 

CSHMBU-1.S7; 

■HO  i 

If  (  (MMX  «•)  AMO  (MMX  <  •) )  TXXX  DO  I 
CBMAMGLB-4 .  17 1 

not 

xr  (mmx  >  •)  nn  oot 

CAMQ-MAMX  */MMX  t 
CAMQ— ATM  (CAM6)  t 
CBMAMGLS— CAMS  t 

not 

If  (MANX  <  •)  Tin  DO  1 
CAM6-MMT*  /MMX  t 
CANO— ATM (CAM)  t 
CSMANCLS-CANGO.llt 

not 

•t 

•t 

00  IA-1  TO  3 1 

If  ( (XX ( IA,  V-#)  MO  (TX(  XA, )  >•)) 

Tin  DO  I 

n 1-3 . 1 4 */2 1 10- XM3H2 ♦ IA-1 1 

nOMI(XQ.C)  -fill 

not 

•t 

xr  ( (XX( XA, ) ■•)  AMD  (TX(XA, )  <  •) ) 


MCI  2 


an  QO) 

IQ-UOB2  +1A-1 ; 

PHI-4. 71) 

STORXX( IQ, 4 )  “PHI ) 

OS; 

•) 

i r  (aim.)  >  it  not  do; 

PHI-Y*<IA,)»/XX<IA,)) 

»hi-«an<phi)  i 

IQ-IH3N2+IA-1 1 
«T0***(IQ,6)«PHI) 

■HO) 

•) 

IF  (X*(IA  ,)  <  •)  THEN  DO) 

PHI-T* ( IA, ) t/XK ( IA, ) ) 
nl-ATAN(PHI)  ) 

PHI-J.14  ♦  HI) 

IQ-IK3M2+IA-1) 

STOK**<IQ,«)-fHX) 

BHD) 

•) 

■HO)  ■ 

*) 

MHO-J.  (!,«,.) I 

skho(x,<)-ipx> 

mou.D-nxtxi 

*«H0(  1,2) -WANT) 

SKHO ( X , 3) -CXHTNOl 
BXH0tX,4) -CZHAM6LE) 

SKH0(X,S) •SIZE) 

■RXHT  IDO; 

•PRINT  MEAHX  HKAHT  CXHTRO  CEHAMGLE  SIZE) 
•) 

•) 

■HO) 

PRZHT  STOREI) 


PROGRAM  LISTING  FOR  THE  MATCHING  PROCESS. 


INPOT:  An  eight  column  data  set  called  GETI.  The  first 
four  columns  correspond  to  the  data  from  the  dictionary 
while  the  next  four  to  the  data  from  the  problem  text. 
The  first  of  the  four  columns  corresponds  to  the 
sequence  number  of  the  feature  Fn,  the  2nd  corresponds 
to  the  ofsight  code  of  the  feature,  and  the  3rd  and  the 
4th  corresponds  to  the  distance  and  angle  components  of 
the  feature  vector. 

A  word  is  picked  and  compared  with  features  in  the 
dictionary  with  the  same  ofsight  code. 

OUTPUT:  Plot 8  for  RVALUE  and  SUMABSDIF  of  the  word 
against  features  in  the  dictionary  with  the  same  ofsight 


PROC  MATRIX; FETCH  OKI;  HX-NRW(Cm)  ; 

WORD-CETI(  127:133,3:8); 

•fruit  word; 

•; 

•> 

•; 

•; 

O-J. (1,3,8) ; 

STORE-J. (1,3,,) ; 

00  M  ID  IB  IT  7; 

I- 1 ;  12-1*1 ; 13-1*2 ; 14-1*3 ; 15-1*4 ; 16-1*5 ; 17-1*6 ; 

0FS1-GETI  ( 1 , 2) -M0KD(1,2) ; 

0PS2-GETI (12,2) -MORS (2, 2) ; 

0PS3 -GET! (13,2) -W0RD(3 ,2) ; 

0FS4-GETX(I4,2) *WORD(4,2) ; 

0FS5-GETKI3 ,2)  -WORD (5, 2) ; 

OFS6-GETI(X6,2)-WORD(6,2) ; 

0FS7-GETI ( 17 ,2) -MORS (7,2) ; 

•j 

•l 

*1 

ir  (I0FS1  n  6)  0>  (0FS2  ME  8)  OR  (0FS3  HE  6)  OR  (0FS4  NE  6) 
OR  (0FS3  NE  8)  OR  (0PS6  ME  8)  OR  (0FS6  ME  •) )  THEM 
00  TO  MEET; 

•; 

•t 

*» 

METRIC  -GETI(I:I7,3)-WORD(  ,3); 

•print  metric; 

METRIC-MS  (METRIC) ; 

•PRINT  METRIC; 

DIF-METRIC (*,) ; 

AMGDIF— GET! (1:17,4) -WORD(  ,4); 

COSDEL— COS ( AMGDIF) iSINDEL-SIN  (AMGDIF) ; 

•t 

IF  (  I  -  128)  THEM  00; 

•PRINT  AMGDIF  COSDEL  SIN  DEL;  END; 

•; 

60M0IFC— COSDEL ( * , ) ;SDMDXFS-SZNOEL ( *, ) ; 
SQAC-SOMDIFC#SOMDIFC*SOMDIFS iSUMDXFS ; 

IF  (1-128)  THEM  DO; 

•PRINT  SOMDIFC  SDMDIFS  SO AC; END; 

SS-SQRT(SQAC) ; 

MM-HORD(  ,2) t 
M— MM(*,) ; 

RVALUE- SSt/M ; 

IF  (1-128)  THEM  DO; 

•PRINT  SI  M;END; 

ST0RE(1,1) -GETI(I4,1) ; 

STORE ( 1,2) -DIF; 

STORE (1,3) -RVALUE; 
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O0//ST0REJ 
NEXT  l 
CSV; 

OUTPUT  0  OOT-O  (RENAME- ( CO  U. -FEATURE  COL2-ABSSUM  00L3-RVM.DE))> 

•; 

•l 

*> 

•PROC  PRIST  DATA-O; 

PROC  GPLOT  DATA-O; 

'PLOT  ABSSOM'PEATDRE-1  RVALOE'PEATDRE-2  /  OVERLAP; 

•PLOT  AB5SUM*rEATURE-l ; 

PLOT  ABSSOM'PEATURE-1/  BZERO  VZERO; 

•PLOT  ABSSUM'RVALOE-PEATURE  /I-H3X; 

SZMBOU  I-HEEDLE  V-DIAMOMD; 

TITLE  ABSOLUTE  DIFFERENCE  AGAINST  FEATURE  NO.  ; 

POOTNOTE1  ONLY  FEATURES  OP  THE  SANE  OFSIGST  INFORMATION  ARE  COMPARED; 
PROC  GPLOT; 

•PLOT  RVALOE*FEATURE; 

PLOT  RVALOB*PEATURE-l/  BZERO  VZERO; 

TITLE  ANGOLA*  CORRELATION  AGAINST  FEATURE  NO.  ; 

FOOTNOTEl  ONLZ  FEATURES  OF  TBE  SAME  OFSIGBT  INFORMATION  ARE  COMPARED; 


DMA.  SUE; 

INPUT  H1M2CDEPGH; 

CAMS; 

•  ••IS  2*3 

3  51*2  S3* 

PROC  G PLOT; 

PLOT  N1*H2-1  C*C-2  E*P-3  G*H«4/0VERLAI; 

symbols  i- jo  m; 

SYMBOL2  t-JOIN; 

SYMBOLS  X-JOIN; 

SYMBOL4  X-JOIN; 

TITLE1  PLOT  OP  A  TWO  DIMENSIONAL  SHAPE  LINE; 

TITLE2  SHOWING  THE  SIZE  VARIABLE  OP  THE  FORM  Ml +M2 -CONST; 

PROC  MATRIX; 

FETCH  SIZE; 

THETA-4; 

OKLTA-C. 2 «•/!••; 

CC-J  . ( IIP , 8 , . )  ; 

DO  1-1  TO  !••; 

CC< 1,1) -3 PTH ETA; 

CC (1,2) —5 PTH ETA; 

CC(I,3) -COS(THETA) ; ; 

CC  ( 1 , 4 )  -  ,5*S IN  (THETA)  ; 

CC(I,C)  -IPS  IN  (THETA)  ; 

CC(I.5)-2«COS(THETA) ; 

CC(1, 7)  — 3PCOS (THETA) ; 

CC(X,S)-1.3PSXN (THETA)  ; 

TH  ETA-TH  BTA+DELTA ; 

END; 

OOTPOT  CC  OOT-CCC  (RENAME- (COL1-N1  COL2-H2  COL3-Z  OOL4-A 
C0L5-B  COL6-C  COL7-D  C0L8-B) ) ; 

DATA  NCC  AND  ;SET  CCC;IP  ( (E  >-  •)  AMD  (  D  >-  •) ) ; 

PROC  GPLOT  DATA- NCC; 

PLOT  M1-H2-1  Z*A«2  B»C-3  D-E-4/0VERLAY ; 

SYMBOLl  I- JOIN; 

SYMBOLS  I-SPLINE; 

SYMBOLS  I-SPLINE; 

SYKBOL4  I-SPLINE; 

TITLE1  PLOT  OP  A  TWO  DIMENSIONAL  SHAPE  LINE; 

TXTLE2  USING  THE  SIZE  VARIABLE  OP  THE  FORM  SQRT(N1  *A*M2  )  -CONSTANT; 
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FAST  FOURIER  TRANSFORM 


and 

RESAMPLING  USING  LINEAR  INTERPOLATION. 

INPUT:  TO  THE  RESAMPLING  USING  LINEAR  INTERPOLATION 
A  data  set  named  SHAPEXY  consisting  of  the  x  and  y 
coordinates  of  the  points  between  which  the  linear 
interpolation  is  to  be  carried  out  before  resampling. 
OUTPUT:  OF  THE  RESAMPLING  USING  LINEAR  INTERPOLATION: 

A  data  set  named  SHAPEl  which  has  three  columns.  The 
first  column  is  the  sequence  number  of  the  sample  while 
the  next  two  are  its  x  and  y  coordinates  respectively. 
The  program  is  designed  to  include  the  input  data  in  the 
resamples.  RES  is  a  temporary  array  in  the  program  which 
stores  the  resampled  values.  Depending  upon  the 
resampling  interval  ISAMP  the  dimensions  of  RES  will 
have  to  be  altered  for  every  run. 

INPUT:  TO  THE  FAST  FOURIER  TRANSFORM  PROGRAM: 

A  data  set  of  the  form  similar  to  SHAPEl  described 
above. 

OUTPUT:  OF  THE  FAST  FOURIER  TRANSFORM  PROGRAM:  Plots  of 
the  real  and  imaginary  parts  of  the  transform, 
normalized  by  the  factor  s  and  a  respectively. 
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WOC  HAT*  IX; 
mCB  SHAPEXT; 

NN— HMM  ( SHAPEXT)  ; 

nm-wm-I; 

CXN-J.  (MM, 2, 8)  ; 

C*H  ( 1 1 NM , )  -SHAPEXT ; 
CXN ( MR, ) -SHAPEXT  ( 1 , )  ; 
DIS-J.  (MM,1,»)  > 

RXS-J. (388,2,8) l 


•CAL  CD  LATE  HI  LINXA*  LXNCTH; 


DU-*; 

do  1-1  to  m  n  li 
j-x+l; 

A-au«x,i)-au(j,i)i 
l-CltM(I,2)-C1tN(J,2)  l 
C-A»A+B*B ; 

C-SQRT(C)  I 
DU-DU+C; 

DIS(I,)-C; 

CKO; 


•PINO  THE  RESAMPLING  INTERVAL; 


ISAMP-OIJ 1/128; 


•FIND  THE  POINTS  IN  BETWEEN; 


KK-t ; 


DO  1-1  TO  NN  BY  1; 
EX-XE+1; 

J-I+l; 

ACL-CHN(J,l)-aiM(I,l)  > 
ACL— ASS (ACL) ; 

C1-C*N( J,l) ; 


PASS 


C2-C*H(J,2)  ; 

xp-aw(i,i>; 

TP-CM(I,2); 

US(CX»)  “CM (I, )  > 

■nT-DXS  ( X, )  »/XSAMP; 

HP-IHT(NINT)  ; 

IP  (ACL  <-  .0001) 

00  E-l  TO  HP  >Z  1; 

n-unt 

xp  (C2  >  i»)  no 

XTT-IO+ISAMP; 

XP  (C2  <  TP)  THU  TTT-TO-ISAKP; 

IXP(IX(1)*XP; 

us(u,2)-m; 

rp-mj 

SO  TO  LAB ELI; 

CUD; 

ELSE 

00; 

H-(CW(J,2)  -CU(I,2) )  #/{CM(J,l)-CU(I,l) )  ; 
C-CM(I,2)-H*OUI(I,l) ; 

DO  M  TO  IP  SI  1; 

KX-KK*1; 

OSQ-XSAHPPXSAHPi 
00Q-  (IM)tl; 

SPA— DBQI/SSQ; 

SPA-SORT(SPA) « 

IP  (XP  <  Cl)  THEN 
XXX-SPA+XS; 

IP  (XP  >  Cl)  THU  XXX-XO-SPA; 

TTT-HIXXX+C; 

XO-XXX; 

IP-TIT; 

KXS ( XX , 1 ) -XXX ; 

CSS (XX, 2) -ITT; 

UO; 

END; 

LAB  ELI : 

END; 


AI-J. (KX,1,1P) j 
RSS(liXX,l)-RSS(l;XX,l)  4-AX; 
AZI-J.(XX,1,1PPI)  ; 
RXS(1:XX,2)  -KES(1>XX,2)  ♦AH; 


KKS-USI1; 
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MSS  1 


OUTWIT  *XS  OOT-SSAPSIXSNAMX-ICOLl-X  COL2-T) ) ; 

MM  3HAPX1  ;  J«_*_;IET  SSAPElXP  (S  HI  •  4  X  HE  •)  * 

MOC  SLOT  MTA-SSAPtl; 

PLOT  Z*l» 

TITLE  PLOT  OP  TIB  SHIR  WOK  PLANS' 

TXTLS2  Or  SXSS-Si 


PSOC  NAT* IX; 

PBTCH  SBAPZ1; 
BN-NaOWtSSAPEl)  » 
TrP-J.(W,2,»)  > 

BX-l+{NN*/2) ; 

Nr-PLOOS(MX) l 
MB-NM-HP; 

PEX-J. <NN.2,S) » 
PET-J.<NH,2.S) i 
KKX-SHAPXK  .2)  | 
IKX-SHAPB1(,3) » 

PSXdsNPf )  —PPT(XEX)  i 
PSZdiNP/ )  ■PPTtXNT)  ; 

•» 

•OPESATXON  POLO; 

•l 

DO  X-l  TO  »  ST  It 
J-NN-Idi 
PSX(J,)-PSX(X.)I 
rsT(j,)»rsT(i,) t 
END; 

•l 

•  OPESATXON  FLIP; 

DO  1-1  TO  M  ST  1; 

J-NM~I+1 » 

PEX(J,2)—  PEX(J,2)  ; 
rSY<J,2)— PST(J,2)  ; 

END* 

*J 

•  OPESATXON  ADO  COMPONENTS; 
TPP(.l)-r*X(d)-PET(,2)  > 
TPP( >2) “FXX( ,2) ♦PET ( >1) l 
AS-J. (1,2,«); 

XXPTS-J. (2,2,4) ; 
rtX-TPP; 

PET-TPP; 

NX-TPPd.i) » 

NT-TPrd.l); 

INS- MU-1; 

P*Xd,)-AS; 

PETd,)-Ai; 

ITISODO-4 ; 

XTXSEVENM ; 

NS-NNI/2; 


M-KS>rLOOK{M«) ; 
ab-tpp<2,I) » 

AI«Tff(2,2)» 

IP  (MS  III  •) 

TW 

00; 

mum,) 

rrzsooo-i; 

BR"TFP(IMS,1)  I 
BI"TFP(IM3,2)  I 
CMS  i 

tr  <m  eq  •) 


DO  I 

ITISEVEM*1 1 
i«-Trr(ms,i)  j 
bi-tpp(imb,2) » 

mi 

QAA»AX#/AI| 

IF  (  At  >-  •)  THEM  QAA-ATAM(QAA) | 

IF  (AM  <  «  AMD  AI  >-  •)  THEM  QAA-3.14I5924  535-ATAIMOAA)  ; 
IF  (AM  <  •  AMO  AI  <  •)  THEM  QAA»3.141392(53S>ATAM(OAA) ; 

OSB-BHt/BIl 

IF  (BB  >-  •)TM8M  QBB-ATAM(QBB) > 

IF  (BB  <  ■  k  BI  >-#)  THEM  QBB-3. 1415 926 S3S-AXAH(QBB) ; 

It  (BB  <  I  I  II  <»)  THEM  OBB>3.141592<5334-ATAM(OBB); 

PHI« (QAA+QBB) 9/2; 

ALPHA- (QAA-QBB) 9/2; 

PHIHT  9RI  AL9BA; 

S-AI9AI+BI9BI; 

J-JQBT( J) ; 

FSIXE-19/S; 


•1MTEH  MULTIPLICATION  BY  EXP  (PHI  ♦  H  ALPHA)  I 


DO  1-1  TO  MX  BY  1; 

IF  (I  -<  HP)  THEM  J-I-li 
IF  (I  >  MF)  THEM  DO) 

IF  (ITIBODD  ■  1)  THEM  J-I-KMi 
IF  (ITISEVEM  >1)  THEM  J-I-SM-1; 
EMO; 

K ALPHA  -J9ALPHA; 

PR IMALIH—PR I ♦RALPH  A ; 

KBX-CQB ( PH IMALPH ) ) 
XEX-«XM{PRXMALPH) ; 
EXPTR(I,1)«KEX| 

EXPTK(1,2) *IEX| 

EXPTt  (2,1)  —IEX; 

EXPTH(2,2)aHEXt 
APP-FHX ( I , ) ; 
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PMC  S 


m(x,)-APPacmiti 

■mi 

ru-ruiriizc; 

•» 

•t 

*» 

•nmt  STANDARD  DfeVXATtOH; 

•l 

•» 

swi*rntmi 

UAV-SAVK,*)  ; 

fflUMSKriM) 
scxzc-sqrt(stsxze) » 
nm  rsxzc  a  uxzzt 
AZQ-SIAPSl(,2i3)  I 
NEAN-AZQ(«-,)*/NNl 

AZQ(  >i)  -AXQ( ,1) -J.  (NN,1)  MEAN  ( ,1)  ; 

AZQ{,2>-AXQ(r2W.  (m,l)MIM(,2)  ; 

SS-AZQMXQ; 

S*A-I*(,+); 

axc-aa*(*,)» 

SCA-NNI/SIG; 

SCA-SQRT(SCA)  , 

std-ii/eca; 
print  SCA  STD  I 
FET-FEXS/STD; 

rue  in  iht  ; 

ootpot  rzx  oor-ru  (rename- <cou-cc  0012-10  j» 

OUTPUT  HI  OCT-rZT ( RENAME- (COLl-CI  C0L2-II) )  ; 

OOTPOT  TFT  OQT-TFF  (RENAME- (C0L1 -ME  CQL2-III) )  ; 

DATA  arCX;J-_H_>SCT  rzx  1 
DATA  KFST;J-_H_;S8T  PET; 

DATA  NTFF;J-.Jt_;SCT  TFF; 

PROC  PRINT  DATA— NFU; 

PROC  PRINT  DATA-HFCT ; 

PROC  PRINT  DATA— MTFF ; 

PROC  PLOT  DATA- NFU  > 

PLOT  RR*J 1 

TZTLC  PLOT  OF  TIC  UAL  PACT  Or  TIC  TRANSFORM  OP  Til  REAL  PART; 

•; 

PROC  PLOT  DATA- NFU; 

PLOT  XR*J ; 

TXTLC  PLOT  OF  TIC  XMMXNARX  PACT  OF  TIC  REAL  PACT  OF  TBE 
TRANSFOCM; 

•; 

PROC  PLOT  OATA-IPET;  ' 

PLOT  CX*J; 

TXTLC  PLOT  OP  TIC  REAL  PACT  OF  TIE  TRANCFORM  OF  TIE  IMAGINARY  PACT; 

•; 

PROC  PLOT  DATA-NFIT ; 

PLOT  XX*J; 

TXTLC  PLOT  OP  III  XNMXNACr  PACT  OP  TIE  TRANSFOCM  or  TIE  XMCXNACT  PACT; 


LISTING  FOR  THE  PROGRAM  USED  TO  ROTATE  AND  SHIFT  THE 

DATA 

INPUT:  Two  column  input  representating  the  data  to  be 
rotated. 

OUTPUT:  Two  column  output  representating  the  rotated 
data 


The  input  and  output  data  set  are  both  called  SHAPEXZ. 


MOC  MATRIX; 
ma  SRAPEXT; 

a-ntarisxAPexx)  > 
srr-i; 

IBIMl 

IXI-1K1; 

CLOOBUS-ttSET; 

TEMX-J.  (H,3,.)  I 

I r  (CLOCWIS-SET)  THEM  00  TO  XXX; 

00  I-l  10  M; 

JXJ-IXI-I; 

TEMPX(I,1)-SHAPEXT(I,1) ; 

TEHPX( JXJ,2 13)  •SRA7CXX ( X<2:3)  ; 

US; 

SIAFBXX*TSMVX; 

XXXt 

FREE  TEtffX; 

OUTPUT  S8AFSXX  OOT-SaAPEXY  <RBIAME-(C0U«J  C0L2-X  C0L3-T)); 


THE  LISTING  OF  THE  PROGRAM  FOR  FINDING  THE  CURVATURE: 


INPUT:  The  input  which  follows  after  the  cards  statement 
consists  of  two  columns  representating  the  x  and  y 
coordinates  of  the  data  points. 

OUTPUT:  The  output  is  a  data  set  named  SXY  consisting  of 
eight  columns.  The  first  column  represents  the  sequence 
number  of  the  data  points.  The  2nd  and  the  3rd  represent 
the  input  data.  The  fourth  represents  the  magnitude  of 
curvature.  The  5th  and  the  6th  represent  the  x  and  y 
components  of  the  Vector  Curvature.  The  7th  and  the  8th 
column  represent  the  location  of  the  critical  points 
which  are  defined  as  the  points  where  the  curvature  is 
two  times  the  mean  curvature. 
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LISTING  OP  THE  PROGRAM  TO  INTERFACE 
THE  ADAPTIVE  LINE  OP  SIGHT  OUTPUT 
TO 

THE  PLOT  ROUTINE. 

INPUT:  Data  from  the  Adaptive  Line  of  Sight  program. 
OUTPUT:  PLOTS  for  showning  the  critical  points  obtained 
using  the  Adaptive  Line  of  sight  method. 
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